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

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

广西灵山地区深部断裂系统的重力资料解释

采用重力资料系统处理,确认广西灵山地区深部存在近EW向、NW向、NE向三组断裂,它们的深度均可达到40 km,近EW、NW向断裂的强烈活动是1936年灵山63/4级地震发生的主要原因。     

基于重力资料的南海及邻区断裂分布特征研究

断裂构造作为重力异常一个非常重要的解释成果,它与能源和矿产资源勘查、构造单元划分和大地构造分区有着密切的联系.有些地质界线以断裂为界,有些能源或矿产资源与断裂伴生.本文利用重力异常,采用归一化总水平导数垂向导数（NVDR-THDR）断裂识别方法推断了南海及邻区主要断裂平面分布位置,首次采用重力异常曲率属性方法反演了南海及邻区主要断裂视深度.在南海及邻区推断断裂57条,其中一级断裂14条,二级断裂43条;断裂走向以NE和NW向为主,其次为NEE、NWW和NNE向以及NNW、EW和近SN向;一级断裂长度大部分为600~1400 km,二级断裂长度大部分为400~900 km;断裂视深度约为2~10 km,其中一级断裂视深度大部分为2~10 km,二级断裂视深度大部分为2~8 km.该研究成果可以为南海及邻区油气勘探、矿产资源勘探以及基础地质研究提供依据.

     

基于category5月球内部结构模型的月球重力固体潮理论值计算

在地球、太阳和其他行星的引潮力的作用下,月球表面上任一点重力值将产生周期性的变化称为月球重力固体潮理论值.研究月球重力固体潮有助于我们更好地了解月球的内部结构、地壳运动以及月球的演化过程,这也为登月任务和月球探测提供了重要的科学依据,帮助我们更好地探索和理解月球的特性和性质.本文阐明了月球重力固体潮形成原因,根据category5月球内部结构模型计算出月球重力潮汐因子,依据Newcomb理论,采用ELP2000-85解析星历和Eckhardt月球天平动模型,提出一个计算月球重力固体潮理论值的方法,选取Apollo11作为测站,以波形图的形式输出该测站一天、一个月、一年的月球重力固体潮理论值,并探索讨论月球重力固体潮与月震的关系.结果表明：该测站一天、一个月、一年的月球重力固体潮理论值,最大值分别为：-1309μGal、-1290μGal、-1277μGal;最小值分别为：-1343μGal、-2349μGal、-2349μGal;与地球重力固体潮(这里选择与Apollo11经度纬度相同的坐标作为测站相比),月球重力固体潮变化相对简单,但波动幅度较大;推测月球的重力固体潮与深源月震之间可...     

基于数据空间和稀疏约束的三维重力和重力梯度数据联合反演

随着重力和重力梯度测量技术的日趋成熟,基于重力和重力梯度数据的反演技术得到了广泛关注.针对反演多解性严重、计算效率低和内存消耗大等难点问题,本文开展了三维重力和重力梯度数据的联合反演研究,该方法结合重力和重力梯度两种数据,将L0范数正则化项加入到目标函数中,并在数据空间下采用改进的共轭梯度算法求解反演最优化问题.同时,本文摒弃了依赖先验信息的深度加权函数,引入了自适应模型积分灵敏度矩阵,用来克服因重力和重力梯度数据核函数随深度增加而衰减引起的趋肤效应问题.为了提高反演计算效率,本文又推导出基于规则网格化的重力和重力梯度快速正演计算方法.模拟试算表明,改进的共轭梯度法可以降低反演的迭代次数,提高反演的收敛速度;自适应模型积分灵敏度矩阵,可以有效解决趋肤效应,提高反演纵向分辨能力;数据空间和改进的共轭梯度算法结合,可以更好地降低反演求解方程的维度,避免存储灵敏度矩阵,有效地降低反演计算时间和内存消耗量.野外实例表明,该算法可以在普通计算机下快速地获得地下密度分布模型,表现出较强的稳定性和适用性.

     

月球和地球正常重力特征差异及成因探讨

月球和地球都非常接近于旋转椭球,基于正常椭球求定的正常重力场占据了月球和地球真实重力场的主要部分,比较月球和地球的正常重力,可以更好地理解这两个天体的重力变化特征.基于Clairaut定理求解月球和地球的正常重力,正常重力随纬度的变化呈现了相反趋势,即地球是两极的正常重力大于赤道的正常重力,而月球是两极的正常重力小于赤...     

利用视线加速度数据恢复月球近区重力异常的研究

本文阐述了利用分步最小二乘方法拟合视线加速度来确定月球重力异常的原理和计算步骤,并利用"Lunar Prospector"扩展任务实测的视线方向加速度数据计算了月球近区的重力异常.最后将计算结果和现今最高阶次的月球重力场模型LP165进行了比较和分析,得出标准差在纬度-75°～75°范围内为±27.354 mGal,75°～85°和-85°～-75°范围内为±61.965 mGal的精度,从而验证了该解算方法的可行性和可靠性.     

应用布格重力异常研究郯庐断裂构造

使用布格重力资料对郯庐断裂带的中段部分（沂沭断裂带）进行了研究. 结果表明, 郯庐断裂带莫霍面及地壳内界面均发生错断,断裂带两侧地壳各界面起伏平缓. 该结果与前人的郯庐断裂带是切穿地壳的深大断裂的认识相一致. 在郯庐断裂带两侧,地壳结构明显不同,西侧沉积层较薄,平均在5 km以下;东侧多数在6 km以上;在断裂带中央沉积物最薄,大约为3～4 km. 断裂带东侧莫霍面埋深浅,大约为33～34 km;西侧莫霍面埋深明显增加,达到36～38 km．反映了莫霍面深度在断裂带附近整体是向西增加的. 郯庐断裂带在重力场分布中则表现为一条宽度较大的线性布格重力异常梯度带.      

博阿断裂西段重复地震识别及断层深部滑移速率研究

本文选取博阿断裂西段周围30 km范围内2010年1月—2015年12月期间9个台站记录到的1272次地震事件的波形记录,利用波形互相关技术从中识别出65组(231次)重复地震,并且基于这些重复地震事件估算该断层的滑移速率.结果显示:该区域重复地震比例为18.2%;博阿断裂阿拉山口—精河段、精河—新源段的平均滑移速率分别约为0.99 mm/a和1.35 mm/a,乌苏—石场段、巴仑台—乌什城段的平均滑移速率分别为1.36 mm/a和2.03 mm/a,这与GPS观测结果基本一致;深部滑移速率的极大值分布在2—10 km深度范围内,其中巴仑台—乌什城段在8 km深度处达到了极大值8.6 mm/a.      

重力和重力梯度数据三维相关成像

本文提出了重力异常三维相关成像方法和重力梯度数据三维相关成像方法,并提出了基于异常分离的三维相关成像方法来提高成像分辨率.通过合成Y型岩脉模型和合成多个直立长方体组合模型的重力异常和重力梯度数据试验分析,验证了本文三维相关成像方法可显示出异常地质体的空间赋存状态和等效剩余质量分布,具有良好的纵向和横向分辨率.     

新近月球重力场模型与地形模型的局部导纳和相关性分析

本文利用三个高阶重力场模型LP150Q、GLGM-3和SGM150j以及嫦娥地形模型CLTM-s01,在频率域内使用固定窗口的方法,研究了不同重力场模型的重力/地形局部导纳谱与局部相关谱的全球分布,以及典型质量瘤盆地重力/地形局部导纳谱与局部相关谱分布的特征.结果表明:加入远月面重力场信息有助于增强重力场模型在中高频段的信号强度;三个重力场模型LP150Q、GLGM-3和SGM150j均适合于作近月面重力/地形局部导纳与局部相关性的分析和近月面地球物理参数的估计;重力场模型LP150Q和GLGM-3较适合作远月面大范围的研究,不适合于作远月面重力/地形局部导纳谱与局部相关谱的分析以及作远月面局部区域地球物理参数的估计,SGM150j较适合于作远月面局部区域地球物理参数的估计;近月面大型质量瘤盆地异常质量的尺度比远月面质量瘤盆地异常质量的大,而深度也比远月面的深.     

The global distribution of gravity wave energy in the lower stratosphere derived from GPS data and gravity wave modelling: Attempt and challenges

Five years of global temperatures retrieved from radio occultations measured by Champ (Challenging Minisatellite Payload) and SAC-C (Satelite de Aplicaciones Cientificas-C) are analyzed for gravity waves (GWs). In order to separate GWs from other atmospheric variations, a high-pass filter was applied on the vertical profile. Resulting temperature fluctuations correspond to vertical wavelengths between 400 m (instrumental resolution) and 10 km (limit of the high-pass filter). The temperature fluctuations can be converted into GW potential energy, but for comparison with parameterization schemes GW momentum flux is required. We therefore used representative values for the vertical and horizontal wavelength to infer GW momentum flux from the GPS measurements. The vertical wavelength value is determined by high-pass filtering, the horizontal wavelength is adopted from a latitude-dependent climatology. The obtained momentum flux distributions agree well, both in global distribution and in absolute values, with simulations using the Warner and McIntyre parameterization (WM) scheme. However, discrepancies are found in the annual cycle. Online simulations, implementing the WM scheme in the mechanistic COMMA-LIM (Cologne Model of the Middle Atmosphere—Leipzig Institute for Meteorology) general circulation model (GCM), do not converge, demonstrating that a good representation of GWs in a GCM requires both a realistic launch distribution and an adequate representation of GW breaking and momentum transfer.     

GOCE引力梯度的频谱分析及滤波

GOCE卫星提供的梯度数据含有非常大的低频误差,如何处理这种误差是GOCE数据处理中最为关键的工作之一.本文根据GOCE卫星的运行情况,首先分析了梯度数据的频率特性,推导了频率与阶次的对应关系;并在此之上,介绍了针对低频误差的滤波方法,即移去恢复和向前向后滤波方法,前者可解决滤波中的低频信号损失问题,后者则主要解决了滤波中的相位漂移问题.最终结果表明:引力梯度的时间频谱与球谐展开中的阶次虽不是一一对应的,但各阶所对应的最大截止频率与阶次却有一定的显式表达.同时也表明,本文所采用的滤波方法是有效的,达到了消除低频误差但保留观测频段信号的目的.     

GOCE引力梯度的频谱分析及滤波

GOCE卫星提供的梯度数据含有非常大的低频误差,如何处理这种误差是GOCE数据处理中最为关键的工作之一.本文根据GOCE卫星的运行情况,首先分析了梯度数据的频率特性,推导了频率与阶次的对应关系;并在此之上,介绍了针对低频误差的滤波方法,即移去恢复和向前向后滤波方法,前者可解决滤波中的低频信号损失问题,后者则主要解决了滤波中的相位漂移问题.最终结果表明:引力梯度的时间频谱与球谐展开中的阶次虽不是一一对应的,但各阶所对应的最大截止频率与阶次却有一定的显式表达.同时也表明,本文所采用的滤波方法是有效的,达到了消除低频误差但保留观测频段信号的目的.     

重震反演中国东北地壳上地幔三维密度结构

本文利用重力和地震P波到时数据反演得到了中国东北地区地壳上地幔三维密度结构.与单一的重力或地震反演相比,重震反演一方面有效地克服了重力反演结果垂向分辨率低的问题,另一方面也提高了地震反演结果的可靠性.结果显示:中国东北地区的地壳及上地幔剩余密度异常分布与构造单元具有明显的相关性,造山带对应低密度异常,盆地对应高密度异常;区域内火山下方有明显的低密度体存在,可能是由于太平洋板块俯冲进入上地幔并部分滞留,在滞留板块深部脱水和软流圈热物质共同作用下产生了上涌岩浆,喷发后形成了火山.     

基于阈值约束的协克里金法联合反演重力与重力梯度数据

常规协克里金方法反演重力或重力梯度数据具有抗噪性好、加入先验信息容易等优点,其反演的地下密度分布能够识别异常体中心位置,还原异常体基本形态,但反演图像光滑,分辨率低,这是由于常规方法估计的密度协方差矩阵全局发散、平稳.为了通过协克里金方法获得聚焦的密度分布需要改善密度协方差矩阵的性质.首先,本文推导了理论密度协方差公式,其性质表明,当理论模型聚焦分布时,其密度协方差矩阵是非平稳且聚焦分布的.为了打破常规协方差矩阵全局平稳、发散的特征,本文设置密度阈值处理协方差矩阵,通过不断更新协方差矩阵来迭代实现协克里金反演,最终得到相对聚焦的反演结果.用本文方法处理重力与重力梯度数据恢复两种密度模型,均得到了与正演模型匹配的反演结果;再将方法运用于文顿盐丘的实际测量重力与重力梯度数据,反演结果与已知的地质情况匹配较好.

     

Main crustal discontinuities of Morocco derived from gravity data

Sharp linear gradients in maps of potential field data are generally assumed to result from sharp discontinuities or boundaries between rocks having different densities or magnetic susceptibilities and are usually associated with faults or other geological contacts. The computation of the horizontal gradients of the gravity field permits us to localize the limits of such blocks and then the fault locations. The horizontal derivative maxima of the Bouguer anomaly and its upward continuation at several heights show lineaments that could reflect the layout of faults and/or contacts and their dip directions. The application of this method to the Bouguer anomaly map of Morocco (with 19,571 points, using an average crustal density ρ = 2.67 g/cm³) allowed us to perform a multiscale analysis of the gravimetric lineaments of the country. The obtained structural map is consistent with several faults already identified in previous studies, and highlights five new major subsurface faults systems with location and dip: the Saghro fault system; Bou-Arfa Midelt fault system; Sidi Slimane Mezquitem fault; Ksar El Kebir–Chefchaouen fault and the Rifan West Mediterranean fault. In addition, this study suggests a new shape and localization for the Agadir-Oujda trans-Moroccan major fault with a NE-SW direction and 900 km length, subdividing Morocco into two main domains. The results of this study contribute to the improvement of the regional structural map of the north western part of Africa, which is situated within the convergence zone between Africa and Eurasia.