青藏高原东北缘地壳地震各向异性研究进展

青藏高原东北缘是青藏块体与华北块体的接触前缘部位,是研究青藏高原隆升扩张和深部动力学问题的重要区域。本文收集了青藏高原东北缘及其邻区由不同方法和不同资料获得的地壳地震各向异性结果,介绍了中上地壳和全地壳各向异性特征;结合区域地质构造、地表运动、构造应力和深部结构,分析了研究区域地壳各向异性的区域分布特征及其与地质构造的关系。结果表明,青藏高原东北缘地震各向异性存在明显的横向区域差异性,体现区域深部构造和地壳介质变形的复杂性;上地壳与全地壳的垂向差异性,反映出该区域可能存在各向异性分层现象。由于青藏高原隆升在其东北缘的伸展边界、物质运移及深部动力模式等尚处在探讨之中,结合多种数据并综合多种方法分析,有助于获得精细、准确的地震各向异性信息,为研究青藏高原隆升演化机制和深部动力模式提供有效的约束。     

Elastic wave velocities in a granitic geothermal reservoir

Geothermal resources have potential for providing cost-effective and sustainable energy. Monitoring of production-induced changes in geothermal reservoirs using seismic waves requires understanding of the elastic properties of the rock and how they change due to injection of fluids and opening and closing of natural and hydraulic fractures. P- and S-wave velocities measured in a granitic geothermal reservoir using sonic logging are systematically lower than those predicted using the composition of the rock. Cracks may occur in granitic rocks from tectonic stresses and from the thermal expansion mismatch between differently oriented anisotropic crystals. An isotropic orientation distribution of microcracks causes a significant reduction in both the P- and S-velocities, consistent with the observed sonic P- and S-velocities. Vertical fractures cause a difference in the velocity of vertically propagating shear waves polarized parallel and perpendicular to the fractures. An assumption that the lower measured velocities are caused by the presence of vertical fractures is inconsistent with the sonic data. This is because vertical fractures cause a decrease in slow S-wave velocity that greatly exceeds the decrease in P-wave velocity, in contrast to the observed data. The growth of vertical fractures in the geothermal reservoir may be monitored using the difference in velocity of the fast and slow shear waves, while the change in P-velocity in a crossplot of measured P- and slow S-velocities is useful for estimating the ratio of the normal-to-shear compliance of the fractures.     

斯通利波和弯曲波联合反演地层VTI各向异性的阵列声波处理方法

利用声波测井获取地层各向异性已成为储层油气评价的一项重要的课题.本文提出了一种单极和偶极振型联合反演地层各向异性的阵列声波处理方法,采用VTI（vertical transversely isotropic）介质斯通利波和弯曲波频散特征,分频段联合求取各向同性面内剪切模量,反演精度相比于传统方法得到了提高;研究了一种同时适用于单极和偶极声波测井的仪器标定方法,消除了仪器对井孔声传播的影响.通过对数值模拟数据的处理和现场测井数据的应用,验证了新方法的正确性和有效性.本文工作为利用阵列声波测井资料评价横观各向同性地层中横波速度的各向异性提供了有效的处理方法.     

Upper crustal anisotropy observed around the Longmenshan fault in the 2013 MS7.0 Lushan earthquake region

Based on the shear wave splitting analysis of the seismic recordings at 17 temporary stations and three permanent stations, we measured the shear wave splitting parameters (i.e., the polarization direction of fast shear wave and the time delay of slow wave) to perform a systematic analysis of the crustal seismic anisotropy around the Longmenshan fault in the 2013 M_S7.0 Lushan earthquake region. We observed apparent spatio-temporal characteristics in the shear wave splitting parameters. The spatial distribution of fast polarization directions showed a clear partitioning in the characteristics from northwest to southeast in the focal region, which changed from NW-SE to NE-SW. In the northwest of the focal region, the fast polarization direction was oriented to NW-SE, which was parallel to the maximum horizontal compressive stress direction. However, the NE-SW fast polarization direction in the southeast of the focal region was parallel to the Longmenshan fault strike. For station BAX on the Central fault in the middle of the focal region, the distribution of fast polarization directions showed a bimodal pattern, with one dominant in the NE-SW direction and the other in the NW-SE direction. With regard to the temporal variation, the time delays were large in the initial stage after the mainshock but then gradually decreased over time and tended to be stable in the later period. This indicated that stress in the focal region increased to a maximum when the main shock occurred, with the stress release caused by the mainshock and aftershock activity, and the stress gradually decreased after a period of time. The scatter of fast polarization directions was large after the main shock, but over time the scatter gradually decreased, indicating that the Lushan earthquake caused a large perturbation in the local stress field. As the stress gradually decreased and was adjusted by the aftershock activity, the perturbation gradually weakened.     

Experimental verification of effective anisotropic crack theories in variable crack aspect ratio medium

Physical modelling of cracked/fractured media using downscaled laboratory experiments has been used with great success as a useful alternative for understanding the effect of anisotropy in the hydrocarbon reservoir characterization and in the crustal and mantle seismology. The main goal of this work was to experimentally verify the predictions of effective elastic parameters in anisotropic cracked media by Hudson and Eshelby–Cheng's effective medium models. For this purpose, we carried out ultrasonic measurements on synthetic anisotropic samples with low crack densities and different aspect ratios. Twelve samples were prepared with two different crack densities, 5% and 8%. Three samples for each crack density presented cracks with only one crack aspect ratio, whereas other three samples for each crack density presented cracks with three different aspect ratios in their composition. It results in samples with aspect ratio values varying from 0.13 to 0.26. All the cracked samples were simulated by penny‐shaped rubber inclusions in a homogeneous isotropic matrix made with epoxy resin. Moreover, an isotropic sample for reference was constructed with epoxy resin only. Regarding velocity predictions performed by the theoretical models, Eshelby–Cheng shows a better fit when compared with the experimental results for samples with single and mix crack aspect ratio (for both crack densities). From velocity values, our comparisons were also performed in terms of the ε, γ, and δ parameters (Thomsen parameters). The results show that Eshelby–Cheng effective medium model fits better with the measurements of ε and γ parameters for crack samples with only one type of crack aspect ratio.     

2017年九寨沟MS7.0地震震源区横波分裂变化特征

对2017年九寨沟M_S7.0地震序列的横波分裂的时空变化特征进行了分析.通过横波窗内S波质点运动图的分析,从九寨沟地震震源区各个地震台站的近震横波记录中提取了横波分裂的快波偏振方向和慢波延迟时间.观测结果显示,震源区各台站的上地壳各向异性在空间上存在分区特征,时间上有随时间的趋势性变化特征.空间上,位于震源区北部余震区内的3个台站中,发震断层东面的台站L5112和L5111只有一个突出的快波偏振优势方向（NNE向）,而西面的L6202台有两个快波偏振优势方向（除了NNE向,还有一个近EW向）,体现了余震区剧烈调整的地壳应力和构造复杂断裂的综合作用;余震区外的3个台站中,震源区东部靠近塔藏断裂（东）附近的JZG台的快波偏振优势方向为NW向,与塔藏断裂（东）的走向一致,南部的台站L5110和L5113的快波偏振优势方向为近EW向,与区域主压应力方向一致;余震区内各台站的平均慢波延迟时间大于余震区外各台站,反映了九寨沟地震孕育过程中余震区的应力积累强于其周边区域.时间上,快波偏振方向在主震后前期离散度较大,随着时间的推移,离散度在后期有逐渐变小的趋势;慢波延迟时间在主震后较大,但随着时间的推移,也表现出逐渐减小,趋于稳定.横波分裂随时间逐渐减小和趋于稳定的变化特征反映了九寨沟地震在孕震中积累的应力,随着主震和余震的发生而导致的应力释放和调整,应力大小和调整幅度逐渐减小,后期趋于稳定.     

Influence of frequency-dependent anisotropy on seismic amplitude-versus-offset signatures for fractured poroelastic rocks

Frequency-dependent amplitude variation with offset offers an effective method for hydrocarbon detections and analysis of fluid flow during production of oil and natural gas within a fractured reservoir. An appropriate representation for the frequency dependency of seismic amplitude variation with offset signatures should incorporate influences of dispersive and attenuating properties of a reservoir and the layered structure for either isotropic or anisotropic dispersion analysis. In this study, we use an equivalent medium permeated with aligned fractures that simulates frequency-dependent anisotropy, which is sensitive to the filled fluid of fractures. The model, where pores and fractures are filled with two different fluids, considers velocity dispersion and attenuation due to mesoscopic wave-induced fluid flow. We have introduced an improved scheme seamlessly linking rock physics modelling and calculations for frequency-dependent reflection coefficients based on the propagator matrix technique. The modelling scheme is performed in the frequency-slowness domain and can properly incorporate effects of both bedded structure of the reservoir and velocity dispersion quantified with frequency-dependent stiffness. Therefore, for a dispersive and attenuated layered model, seismic signatures represent a combined contribution of impedance contrast, layer thickness, anisotropic dispersion of the fractured media and tuning and interference of thin layers, which has been avoided by current conventional methods. Frequency-dependent amplitude variation with offset responses was studied via considering the influences of fracture fills, layer thicknesses and fracture weaknesses for three classes amplitude variation with offset reservoirs. Modelling results show the applicability of the introduced procedure for interpretations of frequency-dependent seismic anomalies associated with both layered structure and velocity dispersion of an equivalent anisotropic medium. The implications indicate that anisotropic velocity dispersion should be incorporated accurately to obtain enhanced amplitude variation with offset interpretations. The presented frequency-dependent amplitude variation with offset modelling procedure offers a useful tool for fracture fluid detections in an anisotropic dispersive reservoir with layered structures.     

Effect of anisotropy on solute transport in degraded fen peat soils

Peat soils are heterogeneous, anisotropic porous media. Compared to mineral soils, there is still limited understanding of physical and solute transport properties of fen peat soils. In this study, we aimed to explore the effect of soil anisotropy on solute transport in degraded fen peat. Undisturbed soil cores, taken in vertical and horizontal direction, were collected from one drained and one restored fen peatland both in a comparable state of soil degradation. Saturated hydraulic conductivity (K _s) and chemical properties of peat were determined for all soil cores. Miscible displacement experiments were conducted under saturated steady state conditions using potassium bromide as a conservative tracer. The results showed that (1) the K _s in vertical direction (K _sv) was significantly higher than that in horizontal direction (K_sh), indicating that K _s of degraded fen peat behaves anisotropically; (2) pronounced preferential flow occurred in vertical direction with a higher immobile water fraction and a higher pore water velocity; (3) the 5% arrival time (a proxy for the strength of preferential flow) was affected by soil anisotropy as well as study site. A strong correlation was found between 5% arrival time and dispersivity, K _s and mobile water fraction; (4) phosphate release was observed from drained peat only. The impact of soil heterogeneity on phosphate leaching was more pronounced than soil anisotropy. The soil core with the strongest preferential flow released the highest amount of phosphate. We conclude that soil anisotropy is crucial in peatland hydrology but additional research is required to fully understand anisotropy effects on solute transport.     

天山构造带地震各向异性及动力学机制研究进展

天山构造带及邻区的深部动力学机制是地球动力学研究的热点,而地震各向异性是区域构造深部动力学机制的一个重要性质。研究表明,天山构造带上地壳各向异性结果呈现区域性分区,受到构造带与断裂走向和区域应力影响;上地幔各向异性的结果认为快波偏振方向和构造带走向基本平行,但在伊塞克湖附近、塔里木盆地和准噶尔盆地挤压区域各向异性快波方向变化复杂,垂直方向上的变化可能由区域性双层各向异性引起,但局部复杂性原因有待进一步探讨。诸多研究支持天山构造带的地壳与上地幔垂直连贯变形机制。此外,地幔柱、软流圈变形、小尺度地幔对流等概念均被用来解释天山构造带的动力学背景,表明该地区的深部动力学机制非常复杂,需要更深入的探讨。     

Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media

A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account.