泥岩、页岩声速各向异性及其影响因素分析

在实验室超声波频率下(纵波主频为700kHz、横波为250kHz)对层理发育的页岩和泥岩的各向异性进行了研究,给出了在干燥和油饱和条件下,样品不同方向上纵、横波速度以及各向异性参数随压力的变化规律. 用X 射线衍射和扫描电镜分析了样品中引起各向异性的原因,指出平行于层理定向排列的粘土矿物和微裂隙是使样品显示出强弹性各向异性的内在原因. 随着压力的增高微裂隙逐渐闭合,样品的各向异性程度减弱. 孔隙流体的存在增强了孔隙（裂隙）的刚度,减弱了各向异性随压力增大而减小的趋势.     

基于面波频散的三维横波速度方位各向异性层析成像方法

地震各向异性是反映地球内部介质特性的重要指针之一。常用的横波分裂法和二维面波方位各向异性层析成像方法很难准确反映各向异性随深度的变化。将与周期相关的区域化面波方位各向异性转换成与深度相关的一维横波速度方位各向异性可以弥补深度信息不足的缺陷。现有三维横波速度各向异性研究多是通过两步方法来实现的,即逐个周期二维面波方位各向异性层析成像以及逐个格点一维横波速度方位各向异性反演。这种分步反演的方式既不利于三维先验约束的引入,也不利于利用原始观测拟合误差对三维模型进行直接评估。因此本文开发了基于面波频散曲线的三维横波速度方位各向异性层析成像方法,并编制了相关正演和反演程序。为了检测方法和程序的有效性,我们对规律分布的三维检测板模型进行了模拟测试。测试结果显示:该方法可以很好地恢复各向同性波速异常、各向异性相对强度和快波方向等三维结构信息;而且反演模型相对于参考模型明显改善了对观测数据的拟合,降低了对观测数据的均方根误差。但对各向同性理论模型进行各向异性反演时,在波速均匀区可产生小于0.5%的假各向异性幅值,在波速非均匀区该假的各向异性幅值会更大,浅部可达3.5%。因此在实际应用中需要谨慎解释（浅部）非均匀区的各向异性结果。      

Investigations of thermal and elastic properties of rocks by means of a cubic press

Investigations of thermal diffusivity and elastic anisotropy have been carried out under high pressure. The pressure is generated by means of a cubic press with a special piston system. The thermal diffusivity for several granulite samples was measured by an impulse method. The elastic anisotropy was investigated for pyrophyllite and calcite. The calcite crystal undergoes a phase transformation in the investigated pressure range with a change of the character of anisotropy.     

Anisotropic 4D seismic response inferred from ultrasonic laboratory measurements: A direct comparison with the isotropic response

Pore-pressure depletion causes changes in the triaxial stress state. Pore-pressure depletion in a flat reservoir, for example, can be reasonably approximated as uniaxial compaction, in which the horizontal effective stress change is smaller than the vertical effective stress. Furthermore, the stress sensitivity of velocities can be angle-dependent. Therefore, time-lapse changes in reservoir elastic anisotropy are expected as a consequence of production, which can complicate the interpretation of the 4D seismic response. The anisotropic 4D seismic response caused by pore-pressure depletion was investigated using existing core velocity measurements. To make a direct comparison between the anisotropic 4D seismic response and the isotropic response based only on vertical velocities, pseudoisotropic elastic properties were utilized, and the two responses were compared in terms of a dynamic rock physics template. A comparison of the dynamic rock physics templates indicates that time-lapse changes in reservoir elastic anisotropy have a noticeable impact on the interpretation of 4D seismic data. Changes in anisotropy as a result of pore-pressure depletion cause a time-lapse amplitude variation with offset response as if there is a reduction in V_P/V_S (i.e., pseudoisotropic V_P/V_S decreases), although the vertical V_P/V_S increases. The impact of time-lapse changes in anisotropy on the amplitude variation with offset gradient was also investigated, and the time-lapse anisotropy was found to enhance changes in the amplitude variation with offset gradient for a given case.     

Seismic Anisotropy in the Deep Mantle, Boundary Layers and the Geometry of Mantle Convection

—An attempt is made to explore the geodynamical significance of seismic anisotropy in the deep mantle on the basis of mineral physics. The mineral physics observations used include the effects of deformation mechanisms on lattice and shape preferred orientation, the effects of pressure on elastic anisotropy and the nature of lattice preferred orientation in deep mantle minerals in dislocation creep regime. Many of these issues are still poorly constrained, but a review of recent results shows that it is possible to interpret deep mantle seismic anisotropy in a unified fashion, based on the solid state processes without invoking partial melting. The key notions are (i) the likely regional variation in the magnitude of anisotropy as deformation mechanisms change from dislocation to diffusion creep (or superplasticity), associated with a change in the stress level and/or grain-size in the convecting mantle with a high Rayleigh number, and (ii) the change in elastic anisotropy with pressure in major mantle minerals, particularly in (Mg, Fe)O. The results provide the following constraints on the style of mantle convection (i) the SH > SV anisotropy in the bottom transition zone and the SV > SH anisotropy in the top lower mantle can be attributed to anisotropy structures (lattice preferred orientation and/or laminated structures) caused by the horizontal flow in this depth range, suggesting the presence of a mid-mantle boundary layer due to (partially) layered convection, (ii) the observed no significant seismic anisotropy in the deep mantle near subduction zones implies that deformation associated with subducting slabs is due mostly to diffusion creep (or superplasticity) and therefore slabs are weak in the deep mantle and hence easily deformed when encountered with resistance forces, and (iii) the SH > SV anisotropy in the cold thick portions of the D" layer is likely to be due to horizontally aligned shape preferred orientation in perovskite plus magnesiowüstite aggregates formed by strong horizontal shear motion in the recent past.     

Insight into craton evolution: Constraints from shear wave splitting in the North China Craton

The multi-episodic tectonic activities from the Precambrian to Cenozoic, including nucleus formation, cratonic amalgamation, and rejuvenation, make the North China Craton (NCC) an ideal natural laboratory for studying craton evolution. Spatial change in the upper deformation records is an important aspect for understanding cratonic formation and rejuvenation. In this study, we performed seismic shear wave splitting analysis using SKS phases from 50 portable stations. Two different methodologies, shear wave splitting measurement and amplitude analysis of transverse/radial components, produced mutually consistent splitting results. These results showed that the seismic anisotropy beneath the Ordos Block can be divided into three subgroups reflecting the tectonic control. Combining these results with those from previous splitting studies in the eastern NCC, we suggest that the Proterozoic amalgamation generated the seismic anisotropy in the boundary zone between the Ordos Block and the Trans-North China Orogen, while the anisotropy in the eastern Trans-North China Orogen and eastern NCC were possibly associated with the lithospheric rejuvenation during the Late Mesozoic to Cenozoic.     

Cosmic-ray vector anisotropy and local characteristics of the interplanetary medium

Variations in the cosmic-ray vector anisotropy observed on Earth are closely connected with the state of the near-Earth interplanetary medium. Hourly characteristics of vector anisotropy for the period 1957–2013, which were obtained by the global survey method from the data of the worldwide network of neutron monitors, make it possible to study the relationship between the cosmic-ray anisotropy and solar wind parameters. In the present work, we have studied the connection between the equatorial component of anisotropy of cosmic rays with a rigidity of 10 GV and the following parameters: velocity and density of the solar wind; density of the interplanetary magnetic field; and cosmic-ray density variations, in which the spatial gradient of cosmic rays in the interplanetary medium is manifested. The characteristics of cosmic-ray anisotropy at various combinations of the interplanetary medium parameters are compared. The possibility of diagnosing the solar wind state from data on the cosmic-ray anisotropy is discussed.     

弹性各向异性介质岩石体波相速度计算

综述弹性波相速度计算的若干严格数值方法和近似解析方法,对数值算例及对不同方法作了比较.结果显示,不同数值方法得到结果一致,其适用于计算任意复杂各向异性介质岩石的相速度,但不便从物理意义的角度去分析物理参数对速度的影响;而近似解析解虽然精确程度各有不同,但是便于从物理意义的角度去进行速度分析,并在此基础上解析求取群速度这...     

SUGGESTIONS FOR A CONSISTENT TERMINOLOGY FOR SEISMIC ANISOTROPY1

Seismic anisotropy is an unfamiliar concept to many geophysicists and the use of misleading and ambiguous terminology has made it more difficult to understand. I suggest here a consistent terminology in which simple expressions have specific meanings similar to their colloquial meanings. It is hoped that use of such language will help to make the increasing number of papers reporting seismic anisotropy more readily comprehensible to the non-specialist. This not a manual of anisotropy, and it is not intended for theoreticians. It is a list of terms which may make anisotropy a little easier to understand for those more familiar with wave propagation in isotropic solids.     

冻土区天然气水合物各向异性储层的AMT响应特征

陆域天然气水合物通常发育于冻土层下方破碎带和岩层裂隙处,其储层会表现明显电性各向异性特征.音频大地电磁法（AMT）能有效探测陆域水合物储藏范围,且不易受高阻冻土层的压制和干扰,可用于陆域天然气水合物探测研究.本文采用AMT对角各向异性二维正演方法,对多种陆域天然气水合物各向异性储层模型及相关参数进行了模拟试算,分析其视电阻率和相位响应特征.结果表明AMT能清晰显示各向异性储层空间位置和分布情况,并对储层的水合物饱和度变化以及储层数量等特征都有所反映.天然气水合物电性各向异性模拟研究为这一新型潜力能源的勘查工作提供了新的思路.     

岩石电阻率图像及各向异性变化的实验研究

在两组人工样品自由表面以中心点为基准对称布设3条辐射状测线,对样品实施单轴应力加载和卸载后,利用电阻率层析成像方法构建了相应的视电阻率相对变化图像,并计算和绘制了表征裂隙产生和发展速率的视电阻率各向异性系数λ*以及表征裂隙产生和发展方位的各向异性主轴方位角α随应力和深度的变化曲线.结果表明:所有测线所对应的RRC图像均随着应力的变化呈现出相同的变化趋势,即在加载阶段,随着应力的增加,视电阻率相对变化图像中电阻率降低区域逐渐收缩,而电阻率升高区域逐渐扩张,在卸载阶段,随着应力的减小,电阻率降低区域继续收缩,电阻率升高区域继续扩张;样品中的高阻体对其所在部位及附近区域的电阻率增幅有较大影响,而对横越高阻体测线的视电阻率相对变化图像的趋势性变化无影响;对于原始电性为各向异性的样品,随着应力的增加,其各向异性程度降低;裂隙主要在岩样的浅部产生和发展,而在较深部位的裂隙产生和发展的速率相对较低.上述结果有助于解释和理解地震、火山活动和大型地质构造运动引起的视电阻率及其各向异性的变化特征,电阻率层析成像方法可能成为目前地震电阻率观测方法的有益补充.      

利用噪声研究地下介质速度结构及其变化在中国大陆的应用

背景噪声层析成像技术已被广泛应用于地壳和上地幔速度结构的研究,这种方法不依靠地震的发生和人工源爆破,只需记录连续的噪声信号而无需产生信号,因为噪声穿过地下介质时会携带信息,然后通过利用台站记录到的连续背景噪声数据进行互相关计算和叠加,即可得到台站间的经验格林函数,从而获取对地下结构的认识。这种方法已经很好地应用于中国的东北地区、华北克拉通、青藏高原以及华南地区,并成功地揭示了这些地区地壳与上地幔顶部的速度结构。此外近年来,一些学者还利用噪声互相关技术研究地下介质地震波速度随时间的变化,通过对汶川大地震前后连续噪声记录的研究发现,大震发生后呈现同震波速降低和震后波速逐渐恢复的特点,这表明可以通过观测地震波特性的变化来监测地下应力的变化,从而为大震的预测预防工作提供科学依据。本文主要综述了近些年来背景噪声技术及其在中国大陆地区的应用。     

Analysis of elastic anisotropy of tight sandstone and the influential factors

Tight sandstone has a certain anisotropy. Using ultrasonic measurements of samples in three different directions and related matched experiments, this study systematically analyzes the pore structure and anisotropy of tight sandstone samples obtained from oil fields and compares results with those of shale. Results firstly show that the anisotropy of tight sandstone is mainly related to the compositional layering and thin interbedding which occur in different sedimentary environments. Tight sandstone has typical transverse isotropic medium characteristics, Young’s modulus increases in different directions with increasing confining pressure, Poisson’s ratio change is not obvious, anisotropic coefficients decrease with increasing effective pressure, and a certain linear relationship exists between ε, γ, and δ. This article finally summarizes anisotropy in different areas, thereby providing a foundation for the use of suitable appraisal models in different regions. This research can be used as an experimental reference for logging evaluation, seismic data interpretation, and fracturing develop of tight sandstones.     

回线源瞬变电磁法有限体积三维任意各向异性正演及分析

目前,瞬变电磁法(TEM)数据基本都是基于各向同性模型进行反演解释,这对于存在明显电性各向异性的勘探区域会产生较大的反演解释误差.为分析电各向异性对回线源瞬变电磁信号的影响方式与程度,本文通过求解离散化的全张量电导率时间域Helmholtz方程,实现了基于有限体积法的TEM任意各向异性的三维正演算法.该算法采用基于交错网格的拟态有限体积法(MFV)对时域Maxwell方程组进行空间域离散,并利用后退欧拉算法(Backward Euler Method)进行时间域离散.为提高时域电磁场的求解精度与效率,该算法将时间分段等步长算法与方程直接求解法相结合.通过对一维各向异性模型以及三维复杂各向同性模型进行测试,验证了本算法对于回线源瞬变电磁响应计算的正确性及有效性.最后,通过对几类典型电各向异性介质中大回线源瞬变电磁信号响应的分析,总结了不同电各向异性类型对TEM电磁信号的影响模式,结果表明,主轴各向异性情况下TEM信号主要受水平方向电导率的影响,倾斜各向异性对TEM信号的影响程度远大于水平各向异性,而通过水平各向异性信号能较清晰判断出各向异性主轴方向.     

SKS波对地壳裂隙各向异性的响应——理论地震图研究

穿透含裂隙、裂缝地壳8s视周期的SV波的理论地震图研究表明,当地壳平均裂隙密度高于0.01即横波各向异性高于1%时,非对称面内不同方位的SKS波均发生分裂;地震图中直接的记录显示是切向T分量上出现SKS波的振动,其振幅随地壳平均裂隙密度的增大而增强,甚至能与径向R分量上的振幅相当.局限于上地壳的强裂缝各向异性同样能引起SKS分裂.长周期SKS波分裂对地壳内裂隙、裂缝的分布缺乏分辨率.直立平行排列裂隙、裂缝使得SKS分裂T分量记录特征具有方位对称性,这来自于HTI介质中快、慢波偏振和到时差随方位变化的对称性;而倾斜裂隙、裂缝使得该方位对称性丧失.对实际观测SKS分裂的偏振解释需要考虑地壳裂隙各向异性,特别是断裂附近的强裂缝各向异性.     

华北克拉通中部过渡带SKS波分裂研究:鄂尔多斯东南角的局部软流圈绕流

本文对布设在华北克拉通东西两块体交界区域的宽频带流动地震观测台阵和部分固定台站的远震波形记录开展了SKS波分裂研究.结果显示,鄂尔多斯块体内部的各向异性比较弱,剪切波分裂导致的时间延迟一般小于0.7s.鄂尔多斯块体东缘的山西断陷带和太行山以及华北平原西部均表现出了比较强的各向异性,时间延迟大于1.0s.特别是在太行山地区观测到的ENE趋向的快波偏振方向明显不同于鄂尔多斯块体和华北平原地区的近E-W和ESE方向的快波偏振方向.在华北克拉通东西两块体交界过渡带的太行山地区观测到的显著上地幔各向异性及变化可能对应于围绕鄂尔多斯块体东南角的局部软流圈绕流,而后者可能起因于鄂尔多斯块体的逆时针旋转以及青藏高原软流圈沿秦岭大别造山带向东的流动.     

THE EFFECTS OF ANISOTROPY ON SQUARE ARRAY RESISTIVITY MEASUREMENTS*

An account is given of the use of the square array technique in investigating the surface effects of rotational anisotropy when the axes of anisotropy are inclined to the surface. It is shown that, as with other arrays, two anisotropy parameters and n can be derived by varying the array orientation. On the basis of these considerations, the effects of such anisotropy on the values of the mean apparent resistivity and azimuthal inhomogeneity ratio normally obtained in square array measurements is reviewed. Particular attention is paid to the variation of resistivity with orientation and it is noted that, in areas of moderate anisotropy, this variation is lower for the square than for the Schlumberger array. In addition to this advantage, the azimuthal inhomogeneity ratio obtained from square array measurements may be used to indicate the severity of anisotropy in an area and two field examples of this use are given. Where anisotropy is severe, gross variations of apparent resistivity with orientation are obtained with either square or collinear arrays. In these circumstances, the use of crossed measurements is considered and the particular stability of the crossed square array demonstrated.     

通过SKS分裂分析非洲中东部 地区地震各向异性

利用非洲台阵(Africa Array)最新的地震记录,通过测量远震SKS震相的分裂参数,详细分析了非洲中东部地区地球介质各向异性可能的成因,包括随应力场变化定向排布的裂隙和岩浆透镜体,以及橄榄石晶格的定向排布等. 结果表明,现今上地幔流动导致的橄榄石晶格定向排布是其各向异性的主要成因,该结果与250 km深度的地幔流动模型一致. 对于少数台站出现的异常结果,则尝试用D″各向异性和双层各向异性模型来解释,并在此基础上讨论了D″各向异性的研究意义.      

上地幔各向异性分析中地壳的影响

应用伪谱法研究了地壳各向异性对上地幔各向异性分析的影响. 用几种不同情况下的两层各向异性模型,模拟了全波场,着重分析了不同深度观测点上的剪切波分裂特征. 结果表明,当地壳为各向同性或其各向异性主轴与上地幔重合或相互垂直时,在地壳观测到的快、慢剪切波的偏振方向保持了上地幔中的偏振方向;重合时,快、慢波走时延迟增大;垂直时,走时延迟减小. 当地壳的各向异性主轴与上地幔成一般夹角时,其粒子运动轨迹变得十分复杂,不能直观地识别上地幔各向异性中的快波方向,对走时延迟也有一定的影响. 将常用的两种分析方法应用于对模拟结果的反演,发现当地壳各向异性与上地幔各向异性方向为一般夹角时,反演结果会出现较大的误差或误判. 进而提出解决这一问题的方法,一是综合多种方法的结果;二是由拟合残差或相关系数等值线作进一步的确定.     

Anisotropy of the upper mantle in Chinese mainland and its vicinity

IntroductionAnisotropy is a universal phenomenon in the upper mantle. The mechanisms, which cause elastic anisotropy of in-situ rocks, include lattice preferred orientation of the minerals composing the rocks and preferred orientation of magma chambers. Lattice preferred orientation (LPO) of olivine is widely believed to be the dominant cause of the upper mantle anisotropy. The observational results of the upper mantle anisotropy can be explained by tectonic process relating to plate motion.…