裂隙型单斜介质中多方位地面三分量记录模拟

针对裂隙型储集层中更具代表性的各向异性介质模型,即在各向同性背景介质中含有两组斜交的垂直裂隙所构成的单斜各向异性介质模型,利用时间和空间上可达任意阶的高阶交错网格有限差分技术,对具有不同裂隙填充物性质的单斜介质中波的传播快照进行了模拟.结果证实各向异性介质中波的传播速度随传播方向的不同而产生明显的差异;裂隙填充物的性质对于速度各向异性具有很大的影响.另外,利用坐标旋转法,对水平层状各向异性介质中多方位地面三分量记录进行了模拟,结果表明了方位各向异性介质中,波的传播速度不仅随入射角的变化而变化,同时也随观测方位的不同而产生差异.数值模拟结果为进一步利用地面多方位地震属性进行各向异性参数的反演及裂隙参数的描述提供理论基础.     

伪谱法地震波正演模拟的多线程并行计算

地震波的正演模拟,尤其是3D正演模拟,往往涉及大规模的数据存储和计算,问题的规模往往超出计算机的物理内存,或者计算时间让问题的求解者难以忍受,即使采用目前存储和计算能力很强的计算机,其计算费用仍然是十分昂贵的.本文提出一种基于多线程协同使用多CPU和计算域分割的正演模拟并行计算技术,使得问题的求解过程得以加快,大大地缩短了用户等待的时间.为了检验我们的并行算法的可行性,文中以傅利叶正演模拟技术为例,给出了声波和3D各向异性弹性波模拟的例子,并对不同版本（串行、并行）运算效率进行了比较,证实了方法的有效性.     

一种弱各向异性介质地震波群速度的近似表示新方法

地震各向异性介质的群速度是关于相角的复杂函数,将其表示成射线角形式较为困难,这给地震各向异性分析以及走时正演模拟等带来诸多不便;另一方面,观测资料表明实际地球介质的地震各向异性通常较弱,这为用射线角近似表示地震波群速度提供了可能.本文基于以射线角近似表示相角的思想,提出了一种弱各向异性条件下,群速度射线角近似表示的新方法.计算表明,在弱地震各向异性条件下,新方法在很宽的射线角范围内,对三种地震波的群速度都能很好地近似,在准SV波计算精度方面显著优于目前通常使用的近似方法.     

Vp/Vs Anisotropy and Implications for Crustal Composition Identification and Earthquake Prediction

Abstract: The ratio of P- to S-wave velocities (Vp/Vs) is regarded as one of the most diagnostic properties of natural rocks. It has been used as a discriminant of composition for the continental crust and provides valuable constraints on its formation and evolution processes. Furthermore, the spatial and temporal changes in Vp/Vs before and after earthquakes are probably the most promising avenue to understanding the source mechanics and possibly predicting earthquakes. Here we calibrate the variations in Vp/Vs in dry, anisotropic crustal rocks and provide a set of basic information for the interpretation of future seismic data from the Wenchuan earthquake Fault zone Scientific Drilling (WFSD) project and other surveys. Vp/Vs is a constant (Ф0) for an isotropic rock. However, most of crustal rocks are anisotropic due to lattice-preferred orientations of anisotropic minerals (e.g., mica, amphibole, plagioclase and pyroxene) and cracks as well as thin compositional layering. The Vp/Vs ratio of an anisotropic rock measured along a selected pair of propagation-vibration directions is an apparent value (Фij) that is significantly different from the value for its isotropic counterpart (Ф0). The usefulness of apparent Vp/Vs ratios as a diagnostic of crustal composition depends largely on rock seismic anisotropy. A 5% of P- and S-wave velocity anisotropy is sufficient to make it impossible to determine the crustal composition using the conventional criteria (Vp/Vs≤1.756 for felsic rocks, 1.7561.944 fluid-filled porous/fractured or partially molten rocks) if the information about the wave propagation-polarization directions with respect to the tectonic framework is unknown. However, the variations in Vp/Vs measured from borehole seismic experiments can be readily interpreted according to the orientations of the ray path and the polarization of the shear waves with respect to the present-day principal stress directions (i.e., the orientation of cracks) and the frozen fabric (i.e., foliation and lineation).     

井筒声场研究回顾与展望

声波测井是测井学科的三大方法之一,在储层评价中起到了重要作用．随着陆上油气勘探目标趋向孔隙结构复杂、非均质性强的油气储集层,基于地层各向同性介质的测井理论显露弊端;因此研究复杂地质条件下的井筒声场分布特征、为储层解释和发展测井新方法提供理论依据成为必然．本文对井筒声场研究进行了回顾,并就当前井筒声场研究现状进行了总结,对其发展趋势和应用前景进行了分析。     

非均质与各向异性黏土基坑抗隆起稳定分析

经典的Prandtl 机构不能很好地反映支护结构入土深度的影响,也不能反映土体强度非均质和各向异性的影响。根据塑性极限分析上限定理,基于修正的Prandtl机构,推导了既可以合理反映支护结构入土深度影响,又可以反映土体强度非均质和各向异性特性的基坑抗隆起稳定的上限解公式,同时基于商业化程序的二次开发,建立了可以考虑土体强度非均质和各向异性影响的强度折减弹塑性有限元方法。通过提出的修正Prandtl机构上限解与强度折减有限元方法以及多块体上限解的对比,验证了修正机构上限解以及强度折减有限元法的合理性和可行性,最后讨论了非均质系数、各向异性系数、支护机构入土深度等因素对基坑抗隆起稳定性的影响,并通过了实例的验证,相关结论可供工程参考。     

大渡河丹巴水电站长大深埋软岩各向异性及宏观力学参数研究

研究了层状岩体各向异性以及宏观力学参数分析方法,提出基于大型岩土分析软件ADINA各向异性材料模块,采用分区分级宏观力学参数模拟分析方法,结合“代表单元体”理论思想,对丹巴水电站二云片岩进行了各向异性模拟和宏观力学参数分析,研究得到丹巴片岩不同片理倾角宏观力学参数以及各向异性特征变化规律。研究结果显示,丹巴二云片岩尺度增加,平行片理向、垂直片理向变形模量 、 收敛幅度不同步, 变化幅度随倾角增大而增大,收敛速率更高,水平向变形模量 趋势较平缓,泊松比的变化规律恰好与之相反;随片理面倾角增大,丹巴二云片岩力学参数趋于极限值规律明显。     

大型地下水封石油洞库渗流场时空演化特征研究

地下水封石油洞库是利用饱水岩体密封性进行石油储存的方式。在地下水封石油洞库的建设中,由于工程体的出现,改变了区域地下水的补给、径流和排泄条件,干扰了原来平衡的地下水渗流场,为保证地下石油洞库的水封效果,需进行裂隙岩体渗透特性及地下水渗流场时空演化研究。以国内首个在建的大型地下水封石油储库为工程背景,结合现场试验数据分析,运用等效连续介质的方法,采用裂隙岩体各向异性渗透张量,建立三维地下水数值模拟模型,预测不同施工进程时地下水位的变化。分析预测结果表明：在地下洞室开挖过程中,无水幕条件下地下水位逐渐降低,主洞室部分区域出现零水头压力区,无法保证水封性;模拟运营期水幕巷道施加定水压力,地下水位上升至设计高度35 m且趋于稳定,可以满足水封效果要求。研究结果对地下石油洞库的水封性评价具有一定的工程意义和理论价值。     

基于离散裂隙网络模型的节理岩体渗透张量及特性分析

节理岩体几何结构非常复杂,研究其渗流特性对于指导含水岩层稳定性分析具有重要价值。应用离散裂隙网络模型DFN方法,基于VC++6.0软件平台,建立了平面渗流分析方法,分析了节理岩体不同几何分布情况下的渗透率张量特征,通过定义渗流定向性系数对岩体渗流的定向性特征进行了定量分析。结果表明：单组节理岩体渗流具有明显的各向异性特征,渗流定向性随着节理角度变化显著;节理随着节理贯通性增加,节理渗透率呈现对数增加趋势;两组节理情况下,各向异性特征随着节理组间夹角变化;两组节理岩体渗流特征研究中,正交分布下,岩体仍存在各向异性,但渗流定向性系数较低;当节理倾角服从正态分布时,随着节理倾角标准差增大,渗透率增加;两组节理夹角不同时,节理渗透主方向倾角随着夹角增大而相应增大,基本沿两组节理夹角方向的角平分线方向。     

Adaptive finite element modeling of marine controlled-source electromagnetic fields in two-dimensional general anisotropic media

In this paper, we extend the scope of numerical simulations of marine controlled-source electromagnetic (CSEM) fields in a particular case of anisotropy (dipping anisotropy) to the general case of anisotropy by using an adaptive finite element approach. In comparison to a dipping anisotropy case, the first order spatial derivatives of the strike-parallel components arise in the partial differential equations for generally anisotropic media, which cause a non-symmetric linear system of equations for finite element modeling. The adaptive finite element method is employed to obtain numerical solutions on a sequence of refined unstructured triangular meshes, which allows for arbitrary model geometries including bathymetry and dipping layers. Numerical results of a 2D anisotropic model show both anisotropy strike and dipping angles have great influence on the marine CSEM responses.