Sensitivity of the elastic anisotropy and seismic reflection amplitude of the Eagle Ford Shale to the presence of kerogen

The Eagle Ford Shale of Central and South Texas is currently of great interest for oil and gas exploration and production. Laboratory studies show that the Eagle Ford Shale is anisotropic, with a correlation between anisotropy and total organic carbon. Organic materials are usually more compliant than other minerals present in organic‐rich shales, and their shapes and distribution are usually anisotropic. This makes organic materials an important source of anisotropy in organic‐rich shales. Neglecting shale anisotropy may lead to incorrect estimates of rock and fluid properties derived from inversion of amplitude versus offset seismic data. Organic materials have a significant effect on the PP and PS reflection amplitudes from the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface. The higher kerogen content of the Lower Eagle Ford compared with that of the Upper Eagle Ford leads to a negative PP reflection amplitude that dims with offset, whereas the PS reflection coefficient increases in magnitude with increasing offset. The PP and PS reflection coefficients at the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface all increase in magnitude with increasing volume fraction of kerogen.     

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.     

多重孔岩石微分等效介质模型及其干燥情形下的解析近似式

经典的微分等效介质(DEM)理论可用于确定多孔介质的弹性性质,但由于缺乏多重孔DEM方程,其估计的多重孔岩石的等效弹性模量依赖于包裹体(即不同孔隙纵横比的孔或缝)的添加顺序.本文首先从Kuster-Toksöz理论出发建立了Zimmermann和Norris两种形式的多重孔DEM方程.Norris形式的多重孔DEM方程预测的等效弹性模量总是位于Hashin-Shtrikman上下限内,而Zimmermann形式的多重孔DEM方程有时会越界.然后,通过使用干燥岩石模量比的解析近似式,对两个相互耦合的Norris形式DEM方程进行解耦得到干燥多重孔岩石的体积和剪切模量解析式.用全DEM方程的数值解对解析近似式的有效性进行了测试,解析公式的计算结果在整个孔隙度分布区间与数值解吻合良好.对实验室测量数据在假设岩石含有双重孔隙的情形下用双重孔DEM解析公式对岩石的弹性模量进行了预测,结果表明,解析式准确地预测了弹性模量随孔隙度的变化.双重孔(即软、硬孔)DEM解析模型可用来反演各孔隙类型的孔隙体积比,它可以通过实验室测量与理论预测之间的平方误差最小反演得到.砂岩样品的反演结果揭示,软孔的孔隙体积百分比与粘土含量没有明显的相关性.     

基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法研究

裂缝储层岩石物理参数的准确获得对地下裂缝预测具有重要意义,而叠前地震反演是获得裂缝岩石物理参数的有效手段.本文从裂缝岩石物理等效模型的构建出发,从测井数据上估测了裂缝岩石物理参数,通过推导含裂缝岩石物理参数的方位各向异性弹性阻抗公式,探讨了基于方位各向异性弹性阻抗的裂缝岩石物理参数地震反演方法.实际工区地震数据应用表明,基于方位各向异性弹性阻抗的裂缝岩石物理参数反演方法合理、可靠,可以降低裂缝岩石物理参数估测的不确定性,为地下裂缝预测提供有力的依据.     

新疆东准噶尔花岗岩类岩石高温高压弹性波速度及其对地壳结构的约束

应用超声波反射-透射法,在最高压力为1.0 GPa（室温）,最高温度为700℃（1.0 GPa）的条件下对新疆东准噶尔地区的卡拉麦里花岗岩带和野马泉岩体的典型花岗岩类岩石（碱长花岗岩、碱性花岗岩、花岗闪长岩、二长花岗岩和石英闪长岩）的纵波速度（V_P）和横波速度（V_S）进行了测量.结果显示,在常温、压力0.4~1.0 GPa条件下,东准噶尔地区花岗岩类岩石的V_P和V_S均随压力呈线性增加,说明在这个压力段岩石中的微裂隙已基本闭合.室温、1.0 GPa时花岗岩类岩石的V_P是5.79~6.84 km·s^-1,V_S是3.26~3.85 km·s^-1.依据压力与V_P及压力与V_S的线性关系,拟合得到常温常压下花岗岩类岩石的纵波和横波压力系数分别是0.1568~0.4078 km/（s·GPa）和0.0722~0.3271 km/（s·GPa）,V_P0和V_S0分别是5.62~6.47 km·s^-1和3.15~3.75 km·s^-1.恒压1.0 GPa、室温到700℃条件下,花岗岩类岩石的V_P和V_S均随温度的升高呈线性降低,温度系数分别为（-3.41~-4.96）×10^-4 km/（s·℃）和（-0.88~-3.22）×10^-4 km/（s·℃）.利用实验获得的花岗岩类岩石的V_P0、V_S0及温度系数和压力系数,结合东准噶尔地区的地热资料,建立了V_P和V_S随深度变化的剖面.将获得的V_P和V_S-深度剖面与该区地球物理探测结果对比,发现东准噶尔地区的碱长花岗岩、碱性花岗岩、二长花岗岩和部分花岗闪长岩的V_P和V_S与该区上地壳速度吻合很好,同时这几种岩石的平均泊松比也与上地壳泊松比一致,因此我们认为这几种类型的岩石是该区上地壳的重要组成部分.另外,石英闪长岩的V_P和V_S均符合中地壳的速度,可能为中地壳中的一种岩石.     

上海中山医院基坑逆作法施工时间效应分析

软土具有流变性,因此,分析基坑开挖时需要考虑时间效应。逆作法施工时,由于出土条件和结构施工条件的限制,导致施工周期延长,因此,基坑施工的时间效应尤为显著。以上海中山医院工程为例,考虑软土流变效应,建立有限元模型对逆作法基坑的施工过程进行模拟,分析了施工期间各工况的围护和土体的内力与变形性状,并与实测结果进行了对比分析。结果表明,考虑土体流变的分析方法能较好的模拟逆作法施工的时间效应,并合理地反映基坑开挖过程对围护以及周边环境的影响。逆作法基坑施工过程复杂,开挖周期相对较长,应合理安排施工流程和施工速度,以达到控制围护结构及土体变形、保护周边环境的目的。     

介观尺度孔隙流体流动作用对纵波传播特征的影响研究——以周期性层状孔隙介质为例

介观尺度孔隙流体流动是地震频段岩石表现出较强速度频散与衰减的主要作用.利用周期性层状孔隙介质模型,基于准静态孔弹性理论给出了模型中孔隙压力、孔隙流体相对运动速度以及固体骨架位移等物理量的数学解析表达式,同时利用Biot理论将其扩展至全频段条件下,克服了传统White模型中介质分界面处流体压力不连续的假设. 在此基础上对准静态与全频段下模型介质中孔隙压力、孔隙流体相对运动速度变化形式及其对弹性波传播特征的影响进行了讨论,为更有效理解介观尺度下流体流动耗散和频散机制提供物理依据.研究结果表明,低频条件下快纵波孔压在介质层内近于定值,慢纵波通过流体扩散改变总孔隙压力, 随频率的增加慢波所形成的流体扩散作用逐渐减弱致使介质中总孔压逐渐接近于快纵波孔压,在较高频率下孔压与应力的二次耦合作用使总孔压超过快纵波孔压.介质中孔隙流体相对运动速度与慢纵波形成的流体相对运动速度变化形式一致;随频率的增加孔隙流体逐渐从排水的弛豫状态过渡到非弛豫状态,其纵波速度-含水饱和度变化形式也从符合孔隙流体均匀分布模式过渡到斑块分布模式,同时介质在不同含水饱和度下的衰减峰值与慢纵波所形成的孔隙流体相对流动速度具有明显的相关性.     

长英质岩石的流变学特征及其影响因素

在总结岩石变形机制与岩石流变学实验进展的基础上,讨论了岩石流变学数据的重复性。虽然高温高压流变学实验积累了大量的数据,但中、上地壳长英质岩石和早期获得的石英集合体的流变实验数据重复性比较差,而近年来发表的石英、长石的流变学实验数据重复性相对较好。虽然利用经验理论模型,根据端元组分可以拟合两相矿物集合体的流变律,但并不能满足定量确定复杂组分和特殊流变性的长英质岩石流变参数的需要。因此,利用长英质岩石流变参数估计大陆地壳流变强度剖面时,即使在相同地温和应变速率条件下,给出的流变曲线、脆-塑性转化带深度也有一定差别,还需要通过大量实验给出更精细的长英质岩石流变学实验数据。根据近年来流变学实验研究的新进展,讨论了在实验室条件下影响长英质岩石流变的各种因素,重点分析了流体、岩石成分、样品粒度和组构对流变的影响。微量结构水对岩石流变有显著的弱化作用,而熔体对流变的影响与熔体含量和分布相关,只有熔体呈薄膜状湿润颗粒边界时,熔体的弱化作用才显著。成分对岩石流变的影响不仅体现在样品的应力指数等流变参数的变化方面,还体现在样品从半脆性变形向塑性变形的转化温度方面。粒度主要影响岩石的变形机制,其中,细粒样品在扩散蠕变域具有应力与粒度线性负相关特性,是理想的应力计,可以用来定量确定韧性剪切带的流变强度;而在位错蠕变域,应力与粒度没有依存关系,这为将实验室条件得出的流变数据外推估计地壳流变提供了重要依据。组构和各向异性是地壳中岩石存在的普遍现象,但关于层状组构对多相矿物组成的岩石流变影响的研究非常少,需要通过新的实验来深入研究。     

软岩三轴卸荷流变力学特性及本构模型研究

以典型软岩－泥质粉砂岩为研究对象,进行了不同应力水平下的恒轴压、分级卸围压室内流变试验。试验结果表明,软岩在卸荷条件下的轴向及侧向流变变形较大,各向异性显著。在恒定围压较高时,轴向流变变形较大,而随着围压逐级降低,侧向流变变形发展较轴向更快,试样在破裂围压下的侧向变形远大于轴向。通过深入分析软岩卸荷流变试验成果可知,线性Burgers流变模型能较好地反映各级围压下流变曲线的衰减流变及稳态流变阶段,但对于破裂围压下的非线性加速流变阶段无法描述。基于上述分析,建立一个非线性损伤流变模型,采用Levenberg-Marquardt非线性优化最小二乘法对试验结果进行拟合,并获得软岩的非线性流变参数。经比较,拟合计算与试验曲线吻合程度很高,说明该流变模型能较好地反映软岩卸荷流变各阶段的曲线特征。     

岩体的加、卸载状态与能量的分配关系

岩体是流变体,在高地应力的长时间作用下其内部非均匀应力随着时间自我松弛,逐渐趋向于静水压力状态。在这样的状态下,虽然其体积仍维持弹性,但丧失了进一步抵抗剪切变形的能力,处于流动的状态,其能量储存具有体积的特征。采用平面压缩模型,研究岩石在不同应力水平长期作用下的变形性状及能量分配关系,给出在快速卸荷条件下发生自持续动力断裂破坏条件,以期揭示深部岩体初始地应力水平（能量在体变及形变上的分配关系）对卸荷条件下岩体运动过程中表现的惯性、黏性及弹塑性变形的影响规律,为建立深部岩体加卸荷本构关系打下基础。