Topography of the Moho and earth crust structure beneath the East Vietnam Sea from 3D inversion of gravity field data

The Moho depth, crustal thickness and fault systems of the East Vietnam Sea (EVS) are determined by 3D interpretation of satellite gravity. The Moho depth is calculated by 3D Parker inversion from residual gravity anomaly that is obtained by removing the gravity effects of seafloor and Pre-Cenozoic sediment basement topographies from the free air anomaly. The 3D inversion solution is constrained by power density spectrum of gravity anomaly and seismic data. The calculated Moho depths in the EVS vary from 30–31 km near the coast to 9 km in the Central Basin. A map of the lithosphere extension factor in the Cenozoic is constructed from Moho and Pre-Cenozoic sediment basement depths. The fault systems constructed by the maximum horizontal gradient approach include NE-SW, NW-SE, and N-S oriented faults. Based on the interpretation results, the EVS is sub-divided into five structural zones which demonstrated the different characteristics of the crustal structure.     

Attempt to identify seismic sources in the eastern Mitidja basin using gravity data and aftershock sequence of the Boumerdes (May 21, 2003; Algeria) earthquake

In order to try to identify the seismogenic sources in the epicentral area, we interpreted data collected from gravity and aftershocks in the eastern part of the Mitidja basin after the occurrence of the 21 May 2003 Boumerdes earthquake (Mw = 6.8). The residual gravity anomaly and the horizontal gradient maps revealed the basement shape and density discontinuities. A seismotectonic model obtained from the aftershocks distribution and gravity data is proposed. This model highlights three active faults: one offshore and two onshore. The offshore reverse fault striking NE-SW, parallel to the coast, is consistent with the USGS focal mechanism of the main event, which is assumed to have the most moment release. The two onshore dipping blind active faults are postulated at crossing angles near the SW tip of the main fault. The interpretation is based mainly on the re-location and distribution of aftershocks, and their focal solutions. It is also supported by the basin structures obtained from the inversion and interpretation of residual gravity anomalies, as well as by additional compiled information such as the pattern of coastal co-seismic uplift. This configuration puts forward the failure mode complexity during the main shock. The topography of the basement obtained from 3D gravity inversion shows that all the onshore located aftershocks occurred in the basement, and the area between the two onshore faults rose as a consequence of their sliding.     

郯庐断裂带中段三维密度结构反演与解释

为深入理解研究区的构造运动规律,进一步探求地壳物质可能的运移模式,本文基于Crust1.0地壳结构模型和最新莫霍面深度模型对自由空气重力异常数据进行联合约束,采用三维重力反演技术得到了郯庐断裂带中段的三维密度结构。反演结果显示,研究区的背景场稳定,断裂构造发育部位表现为重力密度低值,特别是营潍断裂带正处于大范围的重力密度低值区并纵贯整个研究区,由此可推断重力密度低值区可指示断裂构造的存在。      

北京南口—孙河断裂带水热系统特征与成因分析

中低温对流型地热资源在华北地区广泛分布,是一种清洁的替代能源.与活动断裂带相关的水热型地热资源是中低温地热系统的重要组成部分.本文基于高精度重力测量、微动测深及钻孔温度测量等数据,从热源、通道、储层和盖层四个方面探讨了南口—孙河断裂带水热系统特征.低重力异常揭示的燕山期花岗二长岩、闪长岩岩体范围为23.8 km~2和14.3 km~2,放射性测井数据计算得到其生热率均值为3.14μW·m~(-3),侏罗系火山岩生热率均值为1.65μW·m~(-3),隐伏岩体和火山岩均难以构成地热系统的附加热源.重力异常显示南口—孙河断裂带宽度约500~800 m,断裂带切割蓟县系雾迷山组白云岩热储层.钻井温度曲线显示断裂带内水热活动强烈,说明该断裂带是导水、导热的重要通道.断裂带南西侧马池口一带第四系松散层与侏罗系火山岩形成了热储盖层,微动测深显示火山岩最大厚度约1500 m.综上源、通、储、盖四个要素分析,该地热系统为热传导一对流复合型,来自京西北山区的大气降水经远距离径流深循环吸收地层热量后沿南口—孙河断裂上移到达裂隙发育的白云岩地层中形成热水.总之,沿南口—孙河断裂带具备了良好的地热地质条件,可达到规模开采的条件.     

川西高温水热活动区深部热结构的地球物理分析

川西高温水热活动区是地热学研究的典型地区之一,该区温泉、热泉点分布广泛,这些温泉大多沿北西-南东向的金沙江断裂、德格-乡城断裂、甘孜-理塘断裂、鲜水河三条主断裂呈条带状分布.对此高温水热活动区开展地热学分析,是研究青藏高原"东构造结"深部地球动力学过程、开发利用川西高原地热资源的重要基础,具有重要的科学研究意义和实际应用价值.本文针对实测氦同位素所发现的壳幔热结构比,依据重力、航空磁法、地震资料,采用反演莫霍面、居里面深度等方法,对川西高温水热活动区的深部热结构进行了分析.同时对中高温热泉密集出露的巴塘、理塘以及康定水热活动区的典型壳幔热流构与深部背景场进行了对比研究.结果表明,地表氦同位素热流估算的深部热流结构与地球物理异常结果保持较好的对应关系.川西地区热流结构中地壳、地幔的热贡献比较接近：巴塘地区平均Q_c/Q=51.38%、理塘地区平均Q_c/Q=54.39%、康定地区平均Q_c/Q=42.42%.川西地区表现出了"温壳温幔"型的深部热结构,但温地壳背景下形成较高地表热流的原因有待进一步研究.     

Gravimetric study of a potential mineral deposit in the Itapororoca region, Brazil

A semi-detailed gravity survey was carried out over an area of 650 km² localized in the Eo-Neoproterozoic coastal zone of Paraiba State where 548 new gravity stations were added to the existing database. Gravity measurements were made with a LaCoste and Romberg model G meter with a precision of 0.04 mGal. The altitude was determined by barometric levelling with a fixed base achieving a 1.2 m measure of uncertainty, corresponding to an overall accuracy of 0.24 mGal for the Bouguer anomaly. The residual Bouguer map for a 7th degree regional polynomial showed a circumscribed negative anomaly coincident with a localized aero-magnetic anomaly and with hydro-thermally altered outcrops, near the city of Itapororoca. The 3D gravity modelling, constrained by geologic mapping was interpreted as a low density, fractured and/or altered material with a most probable volume of approximately 23 km³, extending to about 8,500 m depth. This result is in accordance with a volcanic body associated with hydrothermal processes accompanied by surface mineralization evidence, which may be of interest to the mining industry.     

云南鲁甸MS6.5地震震源区和周边三维密度结构及重力场变化

2014年8月3日云南鲁甸发生M_S6.5地震,震源机制显示该地震为左旋走滑型,断层破裂面以北西向为主,与昭通断裂带近乎垂直.本文基于三维重力反演技术,采用Crust1.0地壳结构模型和最新Moho面深度模型联合约束,反演得到了该地震震源区及周边的三维密度结构特征.反演结果表明：鲁甸地震震源区南北两侧块体运动受到昭通断裂不同程度阻隔,在地质时间尺度上,昭通断裂北段地壳物质运移的阻挡能力相比南段更强.通过进一步调查距离震中较近的流动重力测点的点值变化发现,从2014年3月至2014年6月大海村（P1）重力点值减小,减小幅度达到19 μGal,而位于昭通断裂东北侧的江底（P5）重力点值同期变化显示为增加,其他距离较远测点在震前重力变化不明显.综上结果对鲁甸地震的成因机制给出了初步解释.本文结果对于认识孕震区深部结构及地球物理场时空变化特征具有一定的意义.对识别与局部重力变化有关的地震前兆也同样具有参考意义.     

珠江三角洲地区重力资料解释的深部断裂系统

对珠江三角洲地区布格重力资料系统处理、研究、解释,确认珠江三角洲地区地壳深部存在近EW向、NW向、NE向三组断裂。近EW向断裂是深部主要断裂,大部分断裂深及30km以上,其中有些断裂是以往很少为人们关注的。NW向断裂分布也比较广泛,但深度相对较浅,大都在10km左右,不超过20km。NE向断裂虽然地表出露最广,规模最大...     

Shallow/upper crustal shear wave structure of the Tehran region (Central Alborz,Iran) from the inversion of Rayleigh wave dispersion measurements

Surface wave dispersion curves from microearthquakes are used to obtain group velocity dispersion maps. The calculation of the local dispersion curves for each grid point from these maps then produces the input data to retrieve the 3D shear wave velocity model of the Tehran region. The group velocity maps indicate that the tomographic results agree well with the three main tectonic features and the geological units in the study area. The tomographic maps generally possess high-velocity structures across most of the mountain belts (Central Alborz and east-southeast mountains), whereas the Tehran Basin correlates to a low-velocity structure. Increasing the period in the study area highlights four independent low-velocity zones that reflect faults and fault junction systems. The shear wave velocity profiles indicate that the depth to bedrock exhibits southward variation ranging from ~?300 m to ~?1500 m. We also focus our analysis on the existence of faults within the shear wave profiles and discuss the low shear wave velocity anomalies deeper than 2 km result from the main fault structures (e.g., North Tehran, North-South Rey and Parchin). Furthermore, we argue that the dip angle of the North Tehran fault varies along fault strike, whereas the North-South Rey fault possesses a constant dip angle. Moreover, initial model uncertainties and checkerboard resolution tests are used to identify reliable and robust anomaly features in the 3D shear wave velocity model and 2D tomographic maps, respectively. Microearthquake analysis provides an effective approach for studying the upper crustal structure heterogeneity, especially the fault structure, of the Tehran region.     

重力剖面和激电测深联合勘测在涡河断裂研究中的应用

为查明涡河断裂的位置、性质等参数,选择重力剖面和激电测深联合勘测方法,对指定区域进行物探勘测工作。结果表明:由密度差异引起的重力异常位置和岩石电阻率、极化率异常位置基本吻合,说明以重力剖面和激电测深为手段的联合勘测是查找断裂的有效方法。     

北京平原区断裂构造重力异常识别研究

区域重力异常蕴含丰富的断裂构造特征信息.采用小波多尺度分解、归一化总水平导数垂向导数(NVDR_THDR)和剖面2.5D重力异常反演方法对北京平原区高精度区域重力资料进行了处理,获得了主要断裂构造平面位置、长度、规模、汇交关系及深浅延伸特征,分析了利用重力异常识别断裂的效果.结果表明:(1)北京平原区NE向断裂平面延展长度大、连续性强,NW向断裂分段性明显,连续性弱,但NW向断裂对NE向断裂有切割改造迹象.近EW或SN向断裂分布较为局限.(2)主控断裂平面和垂向延伸特征差异明显,顺义断裂、孙河断裂、永定河断裂等为盖层断裂;南口断裂、孙河断裂西段、二十里长山断裂和张喜庄断裂为基底断裂;黄庄—高丽营断裂南段、南苑—通县断裂、礼贤断裂、夏垫断裂和皮各庄断裂东段为地壳断裂.(3)利用NVDR_THDR峰值异常带的连续性、幅值、宽度及错切关系可有效识别断裂的平面展布特征,通过小波变换获得的不同深度等效层异常结合2.5D剖面反演可有效研究断裂在基岩内部的深浅延伸情况,但重力异常的垂向分辨能力弱,断裂在新生界内部延伸特征需结合其他地质资料进一步分析.     

Efficient volumetric extraction of most positive/negative curvature and flexure for fracture characterization from 3D seismic data

Most positive/negative curvature and flexure are among the most useful seismic attributes for detecting faults and fractures in the subsurface based on the geometry of seismic reflections. When applied to fracture characterization and modelling of a fractured reservoir, their magnitude and azimuth help quantify both the intensity and orientation of fracturing, respectively. However, previous efforts focus on estimating only the magnitude of both attributes, whereas their associated azimuth is ignored in three‐dimensional (3D) seismic interpretation. This study presents an efficient algorithm for simultaneously evaluating both the magnitude and azimuth of most positive/negative curvature and flexure from 3D seismic data. The approach implemented in this study is analytically more accurate and computationally more efficient compared with the existing approach. The added value of extracting most positive/negative curvature and flexure is demonstrated through the application to a fractured reservoir at Teapot Dome (Wyoming). First, the newly extracted attributes make computer‐aided fault/fracture decomposition possible. This allows interpreters to focus on one particular component for fracture characterization at a time, so that a composite fractured reservoir could be partitioned into different components for detailed analysis. Second, curvature/flexure azimuth allows interpreters to plot fracture histogram and/or rose diagram in an automatic and quantitative manner. Compared with the conventional plotting rose diagram based on manual measurements, automatic plotting is more efficient and offers unbiased insights into fracture systems by illuminating the most likely orientations of natural fractures in fractured reservoirs.     

Changes in magnetic and fractal properties of fractured granites near the Nojima Fault, Japan

Abstract Anisotropy of magnetic susceptibility (AMS) has been used to infer finite strain fabrics in plastically deformed rocks, but there are few studies of magnetic properties in fractured fault rocks. Changes in magnetic and fractal properties of fractured granites from the Disaster Prevention Research Institute, Kyoto University (DPRI) 500 m drilling core towards the Nojima Fault and of the well-foliated fault gouge are described. Fractal analysis of fractured granites shows that the fractal dimension ( D ) increases linearly toward the gouge zone of the fault. In weakly fractured granites ( D = 1.05–1.24), it was found that the degree of AMS correlates positively with the fractal dimension, suggesting a fracture-related magnetic fabric due to fracturing. In strongly fractured granites ( D = 1.25–1.50), weaker, nearly isotropic AMS is found, suggesting erasure by the fragmentation of the magnetic minerals. Within the fault gouge zone, an isotropic AMS fabric was found, as well as twofold increases in magnetic intensity and susceptibility. These changes reflect the production of new magnetite grains, subsequently confirmed by hysteresis studies, which suggests that fault gouge might be regarded as the source of the regional geomagnetic field contrast along active faults. Thus, AMS is clearly a potentially useful tool for inferring the fracturing texture of magnetic minerals in fractured rocks and detecting active faults from the high susceptibility contrast of fault gouge.     

滇西北区断层分数维几何学的研究

本文应用分数维几何学的方法,对滇西北区断裂系进行解析并确立了各断裂系的地位;从该区各典型子区活动断裂与水系的几何学分析和对比,探讨了二者的分布特征及其相关性     

从深部断裂系统讨论肇庆地区的发震构造

对肇庆地区重磁资料系统处理得到的图像进行解释和研究,认为肇庆地区深部存在近东西向、北西向、北东向三组断裂系统,东西向肇庆-广州断裂横贯全区,深及30km以上。北西向肇庆-开平断裂图像反映很清楚,但切割深度相对较浅,在10km左右。北东向吴川-四会断裂规模较大,切割深度在20km左右。其中近东西向是深部主要断裂系统,也是主要活动断裂,它与北西向断裂交汇部位是中强震主要发震构造。     

Coupled Models of Brittle-plastic Deformation and Fluid Flow: Approaches,Methods, and Application to Mesoproterozoic Mineralisation at Mount Isa,Australia

Rock deformation has an important effect on the spatial distribution and temporal evolution of permeability in the Earth’s crust. Hydromechanical coupling is of fundamental significance to natural fluid–rock interaction in porous and fractured hydrothermal systems, and in the assessment and production of hydrocarbon resources and geothermal energy. Shearing and fracturing of rocks can lead to the creation or destruction of permeability when fractures or faults form, or when existing structures are reactivated. Changes in stress orientation or fluid pressure can drive rock failure and create dilating fault zones that have the potential to focus fluid flow, or to breach seals above overpressured fluid compartments. Here, numerical models of deformation and fluid flow related to Mesoproterozoic copper mineralisation at Mount Isa, Australia, are presented that show how changes in deformation geometry in multiply deformed geological architectures relate to changes in dilation patterns, fluid pathways and flow geometry. Coupled numerical simulations of deformation and fluid flow can be useful tools to better understand structural control on fluid flow in hydrothermal mineral systems.     

Full gravity gradient tensors from vertical gravity by cosine transform

We present a method to calculate the full gravity gradient tensors from pre-existing vertical gravity data using the cosine transform technique and discuss the calculated tensor accuracy when the gravity anomalies are contaminated by noise. Gravity gradient tensors computation on 2D infinite horizontal cylinder and 3D ??Y?? type dyke models show that the results computed with the DCT technique are more accurate than the FFT technique regardless if the gravity anomalies are contaminated by noise or not. The DCT precision has increased 2 to 3 times from the standard deviation. In application, the gravity gradient tensors of the Hulin basin calculated by DCT and FFT show that the two results are consistent with each other. However, the DCT results are smoother than results computed with FFT. This shows that the proposed method is less affected by noise and can better reflect the fault distribution.     

龙门山断裂带地壳结构的三维建模

为了实现龙门山区域地质信息的科学管理与共享,推动该区域地球物理资料的统一管理,本文在整理现有二维地球物理探测数据的基础上构造了龙门山断裂带区域的三维地壳结构模型,并结合实测的布格重力异常数据对模型进行了分析。结果表明,地层模型正演的整体重力场与实测布格重力异常基本吻合,初步证明该模型的正确性。局部重力场显示:在该模型的尺度下,沉积层对整体重力场的变化贡献较小;中上地壳有不同程度的隆起和坳陷,与前人研究所揭示的低速异常体吻合;地幔表现为规则的自西北至东南单调上升的重力异常梯度带,与前人反演结果的趋势一致,从而证明了模型的可靠性。最后在重力资料约束下对模型进行了反演,进一步修正了模型。      

山西地震带不同基准的重力场变化

对2011—2017年山西及周边的流动重力观测资料,在绝对点控制下重新做平差计算,从不同尺度分析重力场的动态变化和跨断裂重力段差变化,重点对不同基准的重力场累积变化特征进行比较,结果显示:①文中的平均值基准可以明显削弱单期数据或个别测点对重力场造成的非正常影响;②在山西地震带持续出现了区域性的重力异常,形成了与构造断裂较为一致的重力梯度带;③山西地震带北部测段重力段差的变化幅值最大,南部次之但强于中部,反映了北部的构造活动比较活跃。     

攀西地区重力场特征及地壳密度结构

攀西地区位于峨眉山大火成岩省中西部,构造和岩浆特征显著,地震活动强烈.通过对野外重力测量得到的云县—会东和普洱—七甸两条剖面的高精度重力观测数据进行处理和分析,构建了沿剖面的二维地壳密度结构,其中普洱—七甸剖面与孟连—马龙宽角地震剖面部分位置重合.同时结合区域重力异常特征及下地壳视密度填图结果,得到如下初步认识:红河断裂带是南北地震带南段地区重要的构造分界线,断裂带南北向密度结构和莫霍面分布形态存在较大差异,沿走向构造变化.云县—会东剖面上大姚—会东段下地壳底部存在密度较高的壳幔过渡层,结合研究区下地壳底部壳幔过渡层的密度分布特征,认为该过渡层不是攀西裂谷下的"裂谷垫",而是由岩浆底侵作用造成的.