双三次数值模式与理想Lamb波、龙卷扰动对比模拟分析

通过"双三次曲面拟合—时间步长积分—双三次曲面拟合……",实现成为一种双三次曲面拟合新算法数值模式。它对模式大气作非线性(空间二阶可导)描述,是用数值分析(三次样条/双三次曲面)与计算方法近似求解大气运动非线性偏微分方程组。用Z坐标系高分辨率双三次数值模式,采用原始大气运动方程组和采用欧拉时间积分方案,设置理想扰动,对大气声波Lamb波做数值模拟;采用大气不可压假定,除去声波,对大气重力波龙卷扰动做数值模拟。对比模拟发现:理想Lamb波以声速在水平方向传播,扰动能量频散快,仅能在原地持续数秒钟;理想龙卷扰动则在原地可持续约1 min,且波动内圈(外圈)自下至上始终处于风场辐合(辐散)状态,具冷心结构和气压驻波特征,可因水汽进入扰动与借助凝结潜热释放而得到发展。     

γ测井资料反褶积分层解释法在8411矿床的试验对比

γ测井资料的反褶积分层解释法具有方法简单,精度高,速度快,分辨能力强,用微机进行计算等优点。作者通过用反褶积法对8411矿床γ测井资料的解释,对视特征参数的选择条件,γ强度底数的确定和划分矿层边界的原则,提出了初步看法。将该法的试验结果与传统的平均含量法和矿心分析结果进行了对比,还分析了视特征参数的取值大小对矿层线储量及矿化分层的影响,并讨论了引起线储量相对误差的原因。     

铁矿床深部矿体定位预测的一种新途径 ———以西藏恰功铁矿为例

恰功铁矿是近几年新发现的矿山,矿体的形态、大小、位置、产状、边界等几何特征还没有清楚的认识,为了弄清这些问题,文章首先通过小渡多尺度分析方法将恰功铁矿磁异常资料进行了位场信号分离,然后应用欧拉反褶积方法反演了异常与其场源深度的对应定量关系,最后把位场分离后的异常场源深度叠加起来,就可以对矿山深部找矿提供一种全新的思路.     

甘肃白银厂铜多金属矿田重磁异常及构造特征北大核心CSCD

甘肃白银厂铜多金属硫化物矿田是受陆缘弧环境火山机构及其同生断裂控制的典型火山岩赋矿块状硫化物矿床(VHMS)。为了进一步探讨该矿床重磁场特征及与构造和矿床的关系,通过对区内的重磁场异常数据进行位场分离、小波分析和基于张量数据的三维欧拉反褶积自动确定地质体位置和埋藏深度的定量反演计算。结果表明,研究区布格重力场具有西部高、东部低的特征,其磁场可划分为4个磁场区;局部异常按一定分布规律呈圆形或似圆形正负相间分布,不同尺度的重磁细节异常图在一定程度上反映出引起重磁异常的地质异常体具有一定的延深且分布稳定,而且这些地质异常体具有向深部复合的趋势。经位场分离后正负相间的剩余重磁异常分布范围和分布特征大致反映了白银厂奥陶纪中酸性火山岩建造构造为一个继承性的火山穹窿构造,东、西部不同的重磁异常特征说明东西部的火山喷发具有不同的基底、源区和喷发方式。发育NW向、NE向、近NS向、NEE向等4组断裂构造,近NS向断裂F3、F4、F5、F6、F7形成时间稍晚于NEE向断裂,其与NEE向断裂F1、F2共同形成了研究区棋盘网格状的构造分布特征,这两组断裂为研究区内主要的控岩、控矿断裂。火山机构的分布明显受断裂构造控制,火山口集中分布在深大断裂、大断裂或两组断裂的交汇处。矿田内各矿床及成矿有利地段均处于研究区中部低缓重磁异常场内正负磁异常变化的梯度带内。本次研究为控矿要素研究和开展找矿预测工作提供了丰富、翔实的地球物理资料。     

基于非结构有限元的时间域海洋电磁三维反演

海洋电磁法是一种有效的钻前储层评价手段,可识别出海底构造是否储油气,从而减少干井率降低勘探成本.近年来,频率域三维海洋电磁反演得到了快速发展,但受到空气波的影响,其在浅海环境中无法取得很好的效果.为解决这一问题,本文研究时间域海洋电磁数据三维正反演.正演模拟和计算中,我们选择基于非结构网格的矢量有限元方法.该算法中空间离散采用非结构四面体网格,可很好地拟合复杂海底地形条件和地下复杂结构;而对于时间离散,我们采用无条件稳定的后推欧拉方法,以确保任意时间步长数值计算的稳定性.反演计算中,灵敏度信息采用伴随正演隐式进行计算;同时,依据时间域反演方法的特点采用L-BFGS方法计算模型修正量.通过利用合成数据反演结果证明了本文提出的三维时间域反演方法可用于复杂海底环境,特别是在浅海环境下的有效性.     

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

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

基于一致性理论的Drucker-Prager材料弹黏塑本构模型

系统地论述了一致性理论的主要内容,推导了相关的黏塑性切线模量,并在前人工作的基础上研究了类Hoffman材料的推广条件及其一般性的推广方法;将线性Drucker-Prager模型(简称D-P)推广到一致性的黏塑性模型;给出了推广后模型的向后欧拉算法的有限元离散列式;最后将数值算例和试验成果进行对比,验证了一致性的线性D-P模型的有效性.     

一种表征遥感图像深部层次信息的方法

介绍一种定量地表征遥感地质构造深部层次信息的方法，即利用重磁数据三维欧拉反褶积结果，叠加到遥感图像上，定量地、直接地分辨出遥感构造的深度。并应用到松辽盆地西部斜坡带，效果较好。