K_0状态下含砾无黏性土剪切波速特性试验研究

剪切波速是评估无黏性土力学性质的重要参数。为探究现行剪切波速液化判别方法不适用于砾性土场地的原因,利用自主研发的无黏性土剪切波速与相对密度联合测试系统,以汶川地震液化场地实测资料为参考,通过对福建标准砂和花岗岩圆砾按比例混合配制不同含砾量的砾性土试样在控制相对密度下的K0状态剪切波速试验,研究其含砾量对含砾无黏性土的剪切波速影响规律。试验结果表明:相对密度相同的砾性土剪切波速大于砂土,二者剪切波速比值与砾性土的含砾量成指数函数关系,两者比值最大为1.4;对于砾性土试样,在含砾量40%-60%间存在阀值,当含砾量小于该值,剪切波速随含砾量的增加明显增大,当含砾量大于该值,剪切波速的增长趋势放缓。试验研究指出了现行剪切波速判别方法不适用于砾性土场地的原因,揭示了含砾量与砾性土剪切波速的相互关系,提出了相同相对密度砾性土相对砂土的剪切波速修正公式,为建立新的判别方法提供了理论基础,也为岛礁珊瑚吹填土等宽级配无黏性新工程材料的性能研究拓展了技术手段。     

唐山和巴楚地区液化土动力性能比较研究

对以往曾经发生显著液化的唐山地区和巴楚地区的土动力学性能进行对比研究,包括两个地区砂土相对密度与剪切波速关系、剪切波速与标准贯入击数关系等,以此深入了解两个重要地区土层的动力性能,为研究区域化液化判别方法提供基础。结果表明:现行抗震规范液化判别方法数据来源主要为唐山地区,明显不适用于新疆巴楚地区的液化判别;唐山地区与巴楚地区液化场地砂土取相同相对密度时,两个地区液化场地砂土剪切波速差别很小;但唐山地区与巴楚地区液化层砂土平均标准贯入击数和V_s~N_(63.5)关系曲线差异显著,表明以V_s~N_(63.5)关系区分不同地区砂土抗液化能力是可能的。     

粘性土剪切波速不确定性的统计分析

土体剪切波速具有显著的不确定性。基于全国地震安全性评价工作实测的粘性土的大量剪切波速数据,采用2χ检验方法研究了粘性土不同埋深的剪切波速的概率分布。依据所获得的概率分布,采用相应的统计方法给出了粘性土不同埋深的剪切波速的平均值、最大值、最小值、标准差和变异系数,并给出了95%参考值下限和上限,所获得的结果可用于检验场地剪切波速测试结果的可靠性和粗略估计无实测资料场地土的剪切波速。     

砾性土抗液化能力与其液化判别方法研究进展

从砾性土的定义出发,回顾了砾性土从被认为不会液化到其液化现象引起工程师关注的研究历程及国内外研究进展。总结了含砾量、相对密度等因素对砾性土抗液化能力影响关系的研究现状。介绍了国际上对砾性土场地液化判别的研究成果及各种主流方法的优势与不足;分析了剪切波速方法应用于砾性土场地液化判别的可行性。研究认为:砾性土的液化现象正逐渐被科学界及工程界接受;含砾量对砾性土抗液化能力的影响研究仍然存在较大矛盾,相对密度等影响砂土抗液化能力的因素同样影响砾性土抗液化能力;砂土场地应用的标准贯入、静力触探等液化判别方法不适用于砾性土场地,国际上发展了基于动力触探和贝克尔贯入试验的液化判别方法;剪切波速判别方法在砾性土场地的液化判别中具有优势和潜力,是今后研究的方向之一。     

地下水位对砾性土层剪切波速的影响研究

通过地下水位修正来间接反应有效应力对土层剪切波速的影响的方法被应用于我国剪切波速法砂土液化判别,但现行方法不适用于砾性土层。利用GDS大粒径三轴-剪切波速系统测试不同含砾量,不同相对密度的砾性土在不同应力条件下的剪切波速,研究有效应力对砾性土剪切波速的影响。结果表明:砾性土的剪切波速与有效应力满足幂函数关系,其指数参数稳定,不受土性影响;结合对汶川地震震后调查的砾性土的上覆土层分析,推导出适合砾性土层剪切波速判别法的水位修正公式,给出了简化的线性表达;并分析了地下水位变化对剪切波速测试结果的影响。     

北海市土层剪切波速与埋深的关系

收集了2008年以来北海市城区44个工程场地地震安全性评价中的土层剪切波速测试资料,采用幂函数模型分别对粘土、粗砂、砾砂、其它土类的土层剪切波速与埋深进行加权回归分析,得到不同土类的模型参数及拟合优度指标(判定系数R2和误差标准差σ),并对回归方程和回归系数进行了显著性检验;最后将土层剪切波速预测值与实测值进行对比,验证回归关系式的合理性和适用性。结果表明:北海市城区土层剪切波速与埋深呈现出较好的幂函数关系,拟合优度均在0.81以上;给出的土层剪切波速回归关系式能较好地预测出不同埋深土层剪切波速,特别是在土层中部和底部推测结果准确性更高,其推测结果可供北海市城区缺乏剪切波速测试资料的场地参考使用。     

石家庄市区场地脉动特征及土层剪切波速

根据石家庄市区场地脉动和剪切波速的实例资料,分析了石家庄市区场地脉动的特征以及剪切波速,并给出了场地脉动周期与场地土怪固有周期的统计关系。认为该地区的脉动周期随覆盖层厚度及覆盖层层数的增大而减小,而剪切波速则随深度的增加而增大。     

基于H/V谱比法的乌恰地区场地分类研究

场地条件对于地震动特性影响显著,在工程地震领域中强震记录与之对应的场地条件在研究中十分重要。利用H/V谱比法对新疆乌恰地区12个强震台站场地进行场地分类计算,并与实际钻孔剪切波速的场地分类结果对比分析,H/V谱比法和实际钻孔波速分类结果还是有很好的对应关系,为研究区域内缺乏钻孔剪切波速测试资料的场地进行场地分类提供了参考依据。     

成都地区土层剪切波速与埋深的关系

基于成都市地震安评报告中的场地钻孔剪切波速资料,定性探讨了成都地区Ⅱ类场地常见的10类土的剪切波速与埋深的关系;将拟合优度作为模型的综合评价指标,推荐给出针对不同类型土层的拟合函数模型。基于某工程的钻孔数据,利用推荐的数学回归模型和遗传神经网络模型,与其它研究模型计算得到的剪切波速预测值及实测值进行对比分析。结果表明：...     

航空工业组团阎良片区土层剪切波速与深度的统计关系

基于航空工业组团阎良片区地震小区划项目中大量的剪切波速实测资料,利用最小二乘法,采用常用的幂函数、指数函数、线性函数和一元二次多项式函数4种模型,统计回归了该地区土层剪切波速随深度的变化关系;并以拟合优度为评价标准,通过比对分析和讨论,给出了场地土及其主要岩性剪切波速与深度变化的经验回归公式。结果表明,场地土剪切波速与埋深有显著的相关性,所推荐的经验公式可供无波速资料的场地参考使用。     

孔隙介质中弹性性质及波速摸拟研究评述

本文简短地介绍了近年来提出的一种模拟孔隙岩石弹性性质的方法：网络模型方法。它的计算结果与实验结果很接近。     

裂隙介质的等效弹性模量

详细推导了弹性场的Eshelby椭球包体理论,在此基础上推导了裂隙等效弹性模量的一阶近似式（Hudson公式）。目的有二,一是填补国内献中有关该方面研究的空白;二是阐明有关公式的来龙去脉及物理意义,便于正确的应用。     

Drained-to-undrained transition of bulk modulus in fluid-saturated porous rock induced by dead volume variation

We examine the effect of poroelastic boundary conditions when determining elastic properties of fluid-saturated porous rocks from forced-oscillation laboratory experiments. One undesired yet often unavoidable complication in the estimation of the undrained bulk modulus is due to the presence of the so-called dead volume. It implies that some fluid mass can escape the rock sample under applying a confining pressure perturbation. Thus, the dead volume compromises the undrained state required to unambiguously determine the undrained bulk modulus. In this paper, we model data of recently performed low-frequency (0.1 Hz) measurements. Therein, the dead volume has been systematically varied from 10% to 1000% of the pore volume. For the smallest dead volume, the inferred bulk modulus is close to the Biot–Gassmann undrained bulk modulus. With increasing dead volume, the experimentally inferred bulk modulus approaches the drained bulk modulus. We show that the transition from undrained to drained state as a function of dead volume can be modelled with a 1D poroelastic model for the rock sample-dead volume system with a boundary condition that honours the continuity of the fluid volume flux. We discuss the limitations of the 1D model when applied to data recorded at higher frequencies (up to 100 Hz).     

Topographic elastic least-squares reverse time migration based on vector P- and S-wave equations in the curvilinear coordinates

Elastic least-squares reverse time migration has been applied to multi-component seismic data to obtain high-quality images. However, the final images may suffer from artefacts caused by P- and S-wave crosstalk and severe spurious diffractions caused by complex topographic surface conditions. To suppress these crosstalk artefacts and spurious diffractions, we have developed a topographic separated-wavefield elastic least-squares reverse time migration algorithm. In this method, we apply P- and S-wave separated elastic velocity–stress wave equations in the curvilinear coordinates to derive demigration equations and gradient formulas with respect to P- and S-velocity. For the implementation of topographic separated-wavefield elastic least-squares reverse time migration, the wavefields, gradient directions and step lengths are all calculated in the curvilinear coordinates. Numerical experiments conducted with the two-component data synthetized by a three-topographic-layer with anomalies model and the Canadian Foothills model are considered to verify our method. The results reveal that compared with the conventional method, our method promises imaging results with higher resolution and has a faster residual convergence speed. Finally, we carry out numerical examples on noisy data, imperfect migration velocity and inaccurate surface elevation to analyse its sensitivity to noise, migration velocity and surface elevation error. The results prove that our method is less sensitive to noise compared with the conventional elastic least-squares reverse time migration and needs good migration velocities as other least-squares reverse time migration methods. In addition, when implementing the proposed method, an accurate surface elevation should be obtained by global positioning system to yield high-quality images.     

Characterizing the effect of elastic interactions on the effective elastic properties of porous,cracked rocks

Elastic interactions between pores and cracks reflect how they are organized or spatially distributed in porous rocks. The principle goal of this paper is to understand and characterize the effect of elastic interactions on the effective elastic properties. We perform finite element modelling to quantitatively study how the spatial arrangement of inclusions affects stress distribution and the resulting overall elasticity. It is found that the stress field can be significantly altered by elastic interactions. Compared with a non‐interacting situation, stress shielding considerably stiffens the effective media, while stress amplification appreciably reduces the effective elasticity. We also demonstrate that the T‐matrix approach, which takes into account the ellipsoid distribution of pores or cracks, can successfully characterize the competing effects between stress shielding and stress amplification. Numerical results suggest that, when the concentrations of cracks increase beyond the dilute limit, the single parameter crack density is not sufficient to characterize the contribution of the cracks to the effective elasticity. In order to obtain more reliable and accurate predictions for the effective elastic responses and seismic anisotropies, the spatial distribution of pores and cracks should be included. Additionally, such elastic interaction effects are also dependent on both the pore shapes and the fluid infill.     

对拉科斯特重力仪稳定状态的分析

通过对闰科斯特重力仪稳定状态的分析,认为该仪器从断电至恢复供电后进入稳定状态的时间长短与断电,静置时间长短呈正相关。原因是重力仪核心部分的弹性材料的某些物理性质在不同温度下不尽相同。当工区转换时,仪器不能断电并尽量减少静置时间,以保证测量误差满足要求。     

垂直地震力作用下弹性半空间上扁壳基础的解析解

用布希涅斯克定义的弹性半空间内的垂直位移包括两项积分,除了积分号前面系数的差别之外,第一项积分是单层位势而第二项积分为双层位势。若扁壳基础是正高斯曲率的几何曲面,则壳底与半空间表面间的挤压强度就是半空间表面作用的分布垂直荷载。当越过边界时,双层势位的函数值和单层势位的法向导数值发生跳跃。利用这些性质,本文得出布希涅斯克积分的反演公式,从而避开要求解偏微分—积分方程组的巨大数学困难而易于得出解析解。以椭园抛物面扁壳为例说明本文方法的应用。     

A triple porosity scheme for fluid/solid substitution: theory and experiment

Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects.     

弹性介质各向异性研究沿革、现状与问题

各向异性介质中地震波的传播研究是当今地震学研究领域中的前沿课题之一,同时也是地震学中难题之一.由于地下岩石的各向异性主要表现在:地震波速度随传播方向发生变化;不同类型体波间相互耦合;横渡发生分裂:面波速度频散依赖于传播方向等.薄互层与裂隙定向分布等产生视各向异性,放在石油地震勘探、地震预测和岩石层物理与动力学研究中有极大潜力和应用前景,并受到广泛重视.为此,本文较为详尽地讨论了弹性介质中地震波各向异性研究的沿革,简述了国内外现今取得的主要研究成果以及目前尚存在的和有待解决的一些主要问题.最后对我国在各向异性介质中地震波动、检测和应用等研究提出了一些看法.     

岩石中裂纹对弹性波速度的影响

本文在常温常压下,对岩石中的层理,裂纹或裂缝引起的弹性波速变化特征进行了实验研究,并通过简化的人工裂缝模型,研究了裂缝密度和相对位置的变化对波速的影响取得了一些有意义的结果。