谱矩方法在磁源体深度反演中的应用研究

谱矩方法可以对数据的表面形貌做较为细致的描述.它以随机过程为理论基础,用各阶谱矩及统计不变量等具体的参数表征表面的几何形态,算术平均顶点曲率是一种基于四阶谱矩的统计不变量.通常,埋深不同的场源所引起的磁异常尺度不同,从曲率的角度来理解即为磁异常曲面的弯曲程度不同.因此,本文应用算术平均顶点曲率提取磁异常的几何信息,并将所提取的信息用于场源深度的反演.理论上推导了基于谱矩的球状磁源体和板状磁源体的反演公式,得到了场源深度与磁异常、曲率之间的关系式.结合理论模型计算验证了方法的有效性,并与欧拉反褶积方法进行对比.与传统的方法相比,该方法快速简单,无需调节参数,且有较好的反演精度.最后,将该方法用于塔里木盆地航磁异常的反演和解释中,反演出的磁源体深度可满足区域磁异常数据分析和解释的要求,为克拉通沉积盆地磁异常源的深度划分提供丰富的信息.     

基于剩余异常相关成像的重磁物性反演方法

将场源区剖分成长方体单元,通过采集的重磁数据反演出这些单元的密度或者磁化率变化,勾画出场源的分布图像,这种方式是重磁三维反演的重要方向.重磁相关成像通过计算测量的重磁异常与地下各点在测区上的重磁异常的归一化相关,显示出异常地质体的空间赋存状态和等效剩余重磁物性.该方法计算速度快,方法简单、稳定,但是反演的结果只是在-1到+1之间的等效物性,不能够直接反演剩余密度或者磁化率,并且无法引入已知的地质约束.本文通过对物性模型的正演和实测结果的残差进行相关成像,迭代更新物性模型实现对物性参数的反演过程.模型实验证明该方法相对相关成像不仅能提高分辨率,还能够得到真正的物性参数.     

基于二阶导数的磁源边界与顶部深度快速反演

为实现磁异常反演,首先提出了磁异常非参数快速反演的概念,即无需提供先验信息.在Nabighian提出的磁场通用梯度公式的基础上,推导实现了基于磁异常垂向二阶导数的非参数快速反演(V2D_depth).它不仅可以获取场源的边界信息,同时可以反演场源的埋深.通过与Tilt_depth方法对比,本文方法计算的场源边界更清晰,反演的深度也更接近真实深度,同时较大程度上克服了Tilt_depth方法受叠加异常的影响.理论模型验证了方法的有效性,并通过准噶尔盆地某区块磁异常数据的处理,提取了受强大的区域背景场掩盖的石炭系火成岩产生的弱异常,突出了构造分区、断裂分布等信息,获取了火成岩的位置及埋深参数,为该区的火成岩油气藏勘探提供了有效的处理途径.     

基于交叉梯度约束的重力、磁法和大地电磁三维联合反演

为了降低单一地球物理方法反演的多解性及受噪声的影响程度,本文围绕重力、磁法和大地电磁法开展了三维联合反演的研究.重、磁采用基于对数障碍法的正则化反演算法,大地电磁使用limited-memory BroydenFletcher-Goldfarb-Shanno(L-BFGS)反演算法,引入交叉梯度函数实现了三种物性结构的相互耦合,最终开发出一套重磁电三维联合反演算法,并实现MPI并行加速计算.通过理论模型算例验证了算法的准确性,结果表明:不论是单棱柱体模型还是组合棱柱体模型,联合反演结果相较单独反演对于异常体的空间形态刻画以及物性数值恢复具有较好的提升;单棱柱体模型算例使得异常体的物性参数(密度、磁化率和电阻率)更加接近于真实的物性参数;组合棱柱体模型的联合反演结果不仅仅消除了围岩物性参数的假异常,而且还增强了异常体边界结构的恢复程度.     

基于概率成像技术的低纬度磁异常化极方法

化极转换是磁异常解释的重要基础,为了克服在地磁纬度较低的地区尤其是磁赤道附近化极不稳定的问题,出现了多种化极方法．本文基于概率成像技术提出了一种等效物性的反演方法,实现对地下等效场源的反演成像,取得了对低纬度磁异常稳定化极的效果．化极反演中逐次对剩余异常进行反演成像,实现由概率模型到物性模型的复杂映射,避免了类似反演中需要大型方程组求解等问题,并将概率模型的构制、物性参数的反演和反演评价有机地集成到一起,加速了反演成像的进程,使反演成像速度与目前概率成像的计算速度达到可相比拟的程度．对理论模型和实际资料的计算表明,该方法对低纬度磁异常化极处理是稳定有效的,而且可以较好地压制噪声干扰,能够在噪声干扰条件下进行反演化极．     

基于磁异常的边界特征增强方法对比研究

基于磁异常进行场源边界识别(如断裂划分)是磁法勘探解释工作的重要内容.由于受到斜磁化、场源埋深、异常叠加以及噪声等多种因素的影响,基于磁异常直接进行场源的边界识别往往难以准确识别出场源体的边界,所以,人们研究了各种场源边界特征增强的方法技术,以便于进行场源边界的识别工作.近年来,国内外出现了很多新的边界特征增强方法,但有些方法的边界特征增强效果与已有方法相近,有的方法效果有待深化检验,存在验证方法的理论模型过于简单,未全面细致地分析出它们的优缺点等问题.针对上述情况,本文运用多种复杂、贴近实际的理论模型,深入对比分析了包括解析信号模法和倾斜角法等在内的15种具有代表性的边界增强方法,分析了它们在不同磁化方向、不同场源埋深、不同场源形态、异常叠加和噪声干扰等多种因素影响下的边界增强效果,重点归纳总结出它们的优缺点和适用条件,以为实际应用时的方法选择提供参考.通过模型实验,我们认为解析信号模垂向导数法、总水平导数法、解析信号模倾斜角法是在多种因素影响下,适用性较强、应用效果较好的三种方法.     

重磁异常反演的拟BP神经网络方法及其应用

把神经网络与重磁异常反演理论相结合,提出了用于重磁反演的一种拟BP神经网络方法．基于3层神经网络结构,把隐含层神经元设定为三维空间物性（磁化强度或密度）单元．对实测与理论重磁异常经S型函数变换,采用自动修改物性单元物性值的拟BP算法,反演三维空间的物性分布．利用该网络对理论模型数据和内蒙古某花岗岩体上的航磁资料进行了反演计算,取得了满意的反演效果．     

基于方向tilt-Euler的三维磁数据快速反演

三维磁异常快速反演是磁法勘探进行地质解释的一项重要内容.本文在方向tilt梯度基础上进行导数处理,并进行合理改进,得到了6个改进的方向tilt梯度导数,丰富了不同方向上的磁异常信息.对常规欧拉反褶积方程进行三个方向求导,所构成的方程组与6个改进的方向tilt梯度导数进行结合,从而构建出了两个基于方向tilt梯度导数的快速反演方程组(方向tilt-Euler法),用于反演场源的位置参数.另外,将三个方向导数下的欧拉反褶积进行平方和及均方根处理,得到了场源构造指数的求解公式.单一模型试验表明,6个改进的方向tilt梯度导数异常基本不受磁化方向影响,可从不同方向上反映磁异常信息,方向tilt-Euler能够准确地反演出场源的位置及构造指数.叠加模型试验证实了方向tilt-Euler法相对于tilt-Euler法,具有反演解聚集度更强、连续性更好、反演准确度更高、受噪声影响更小及无虚假反演解等优点.将本文方法应用于内蒙古塔木素航磁资料处理中,获得了丰富的地下磁源分布信息,为了解该地区隐伏断裂构造及岩体分布等提供了有效依据.     

数据空间磁异常模量三维反演

强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果.     

多因复成矿区磁异常偏位解释法

分析了多因复成矿床磁性体的典型特征,总结了该类型矿区磁异常反演解释存在的主要问题,通过提取磁异常基本特征、引入解释偏角和剖面位移等概念,提出了以磁异常共轭解释法为基础的偏位解释法,模型试算及实例分析表明该方法能获得较高的反演精度和符合实际的反演结果,能实现复杂磁异常的反演解释.对比分析认为该方法相对传统方法在一些特殊情况下具有明显优势,是磁法勘探领域中一种新的反演方法,能较好解决多生产实际中一些以往不能解决的重要问题.     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Earthquake Engineering and Engineering Vibration Aims And Scope

正This journal is established by the Institute of Engineering Mechanics(IEM),China Earthquake Administration,to promote scientific exchange between Chinese and foreign scientists and engineers so as to improve the theory and practice of earthquake hazards mitigation,preparedness,and recovery.To accomplish this purpose,the journal aims to attract a balanced number of papers between Chinese and     

INTERNATIONAL RESEARCH AND TRAINING CENTER ON EROSION AND SEDIMENTATION(IRTCES)

正Director:Shangfu Kuang,China Vice-directors:Chunhong Hu,China Duihu Ning,China Guangquan Liu,China The International Research and Training Center on Erosion and Sedimentation(IRTCES)was jointly set up by the Government of China and UNESCO on July 21,1984.It aims at the promotion of international exchange of knowledge and cooperation in the studies of erosion and     

Calendar of events

正~~     

Information for authors

正SCIENCE CHINA Earth Sciences,an academic journal cosponsored by the Chinese Academy of Sciences and the National Natural Science Foundation of China,and published by Science China Press and Springer,is committed to publishing high-quality,original results in both basic and applied research.