重力梯度法的改进及其确定地下构造体参数的研究

被引量 : 0次 | 上传用户:xiao678
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
笔者自2007年开始接触重力梯度法,并对其做了较深入研究。开始,为了减少重力梯度法的求解过程产生的多解性,首先利用重力位异常的垂直二阶导数和水平导数,在大区域内确定断裂带的位置,为重力梯度法创造一个合适的计算区间。然后利用重力梯度法在小区域内去计算其内的地下构造体。虽然这种联合求解法在一定程度上抑制了错误的结果,但对于精确求解地下构造体的位置效果不十分理想;由于没有考虑重力异常的求导过程中产生的高频放大误差,对解释结果的准确性进一步降低;同时该方法主要应用于小区域倾斜台阶的计算上,对直立台阶和大区域上的计算有一定的缺陷。但当时改进后的方法能够有效地剔除一些错误的、模糊的信息,为实际资料的处理解释提供了更可靠的基础数据,能够在未知区数据处理中做出较好的推断,为该地区后续其它方法的解释奠定了基础。为了解决重力梯度法计算中出现的上述问题,笔者博士期间又对该方法作了进一步改进,主要工作有如下两点。一、引入别列兹金提出的圆滑滤波因子常规求导为了削弱误差影响,一般采用对求导结果进行圆滑处理,或者在计算前采用滤波等措施减少误差。但当地下存在其它地质体时,就无法利用以上方法来很好地削弱影响,这会对解得的Vzz和Vzx造成较大误差,致使其参量图变形。本文利用求导过程中添加波数域滤波因子[6],同时选择相关参数的方法,可以在构建空间参量图时,大大削弱各种地质体之间的影响。另外引入滤波因子也可以增加计算的稳定性,且在计算机编程中可使得参数选择更加方便灵活。在实际处理中合理[7]地选择别列兹金滤波因子的指数n和谐波数N,可以把干扰源造成的影响减小。重力梯度法的可靠性与所使用资料的精度有很大关系,比例尺过小,会造成很多细小的地层划分不出来,所以在划分地层时,应尽量选取较大比例尺的数据。实际资料处理时,还应根据比例尺的大小,合理地选择图切剖面上的点距。二、引入状态判别因子由于重力梯度法对精度有很高的要求,在数据质量不高或者地形复杂的地区,它对角点位置的判定结果中就会出现多解性,成为一个病态问题。以往解决病态问题的方法,通常是利用磁、电、震等资料来对重力解释结果进行约束,但这种约束方法往往存在以下三个问题:一是不能及时进行信息反馈,使约束效果降低;二是各种约束条件的引入使逻辑判断更复杂,导致解释过程更加混乱;三是大量资料的应用,造成工作量的剧增。为了解决这种弊端,更准确地确定角点位置,本文引入了状态判别因子对确定角点位置进行优化。计算结果表明,这种引入比常规的确定的角点位置更为精确。同时该方法能在一定程度上反证结论,即:利用角点位置反推状态判别因子f ’,如果f ’值的范围与原f在同一区间,且两者之差不超过f值的5%,则证明结果正确,反之则说明求解过程中存在错误。区别于以往对比电法和磁法资料的验证方法,引入状态判别因子后,该方法不需要在一条剖面解释结束后在进行整体验证,如果解释过程中一旦出现问题,可以随时修改,缩小了计算量,减少了逻辑错误。实践证明在重力梯度法中引入状态判别因子f后,计算精度有了很大提高,并且在计算过程中就可以结合反馈的信息(倾角和深度)进行约束,容易剔除错误的解。同时根据相关震、磁、电资料,也可对地下矿体的位置和产状进行判定。本文在进行了理论证明后,把改进后的方法应用于内蒙古道朗和都格地区[8]。并在研究区内做了各种相关的重力数据处理,把得到的结果与重力梯度法得到的结果做佐证和整合,得到如下认识:(1)利用研究区的水平导数和垂直导数,划分出大区域内地下构造体的边界和断裂,为重力梯度法在小区域内的计算提供前提。同时将导数得到的结果与重力梯度法得到的结果做验证对比。如果相似程度高,则说明结果正确,如果相似程度很低,或者出现相反结果,则需对出现错误的区域进行重新计算。一般出现这种问题的情况有两种,一种为错误区域的地下构造体是由多种地质体组合成的,如果是这种情况,则可以认为改进后的重力梯度法的解释结果为正确的;另一种情况是错误区域是简单构造体,但由于附近有凹陷或者溶洞等存在,使构造体出现多种解释结果,如果遇到这种情况,则需参考其他方法的结果,综合处理。(2)结合欧拉反褶积[9]和视密度反演以及小波多分辨分析方法[10],推断研究区内矿体的位置和分布情况。视密度反演可以大致圈定矿体的分布情况,它与对比重力梯度法得到的地下构造体位置具有很好的对应关系;重力梯度法得到的断裂和构造边界,能有效的补充欧拉反褶积反演出的结果;小波多分辨分析,可以有效地分离出重力异常区域,而重力异常区域即为地下构造体的位置,因此,小波多分辨分析得到的结果,可以与重力梯度法得到结果互相验证。(3)由于研究区内主要是金属矿体,因此本文还收集了电法数据和钻井数据,结合钻井和电法解释结果,可以大致推断出研究区内的主要矿体为黄铜矿、黄铁矿以及伴生有小量的钨矿。而通过重力梯度法得到的地下构造体三维解释图,可以看到,研究区中部高,北部和南部低;有一条北东方向的地下地质体贯穿整个研究区,把研究区分成三部分:西北部、中部和东南部;西北部与大地质体相交部分有一隆起带,结合此地区的电法资料,可以判断该区是一个金属矿区;东南部有一条z字型的小地质体带,同样可以得到这个z个型地质体就是该区的矿体主要集中地。(4)为了验证重力梯度法的精确与否,在研究区内选取了4条剖面,并做了数据处理,把得到的结果与该区的综合解释结果做对比,同时也和没改进的重力梯度法得到的结果进行对比,实验结果表明,在引入状态判别因子和滤波因子后,能够更有效的判断出地下构造体的物性参数。重力梯度法作为一个对延展型的截断板和局部型的断裂等地质构造具有较好解释效果的方法,并且它不需要事先估计和给定地下地质体的密度,因此在理论上和实践中都具有非常重要的意义。在理论上,该方法可以很好地补充现有对断裂等构造体的解释方法,为重力勘探理论解释模型提供更多的选择;在实践上,能够提高现有重力勘探对地下构造体的解释精度。与地震等花费巨大的地球物理勘探方法相比,重力梯度法成本低廉,因此将具有很好的理论研究意义和巨大的社会经济效益。
其他文献
本文从分析民意入手,围绕民意特点探究了公正定罪的实现问题。对于该问题的思考源于最近几年民意对刑事审判影响的加重。民众通过各种途径表达出自己对刑事案件的事实认识和法
左右手复合传输线(cRLH TL)是一种特殊的传输线,它同时具有左手传输线和右手传输线的特征。在较低的频率区间内,这种传输线呈现出左手特性,即相位的传播方向与能量的传播方向相
以山桐子嫩枝为扦插材料,选用当年生半木质化发育良好的枝条,采枝条中上部剪取为插穗,分别用单一激素ABT、NAA、IBA及混合激素配比IBA∶NAA(1∶1)4种处理方式,研究不同激素处理
针对雷达多功能小型化的要求,提出一种基于软件化的雷达设计方法,并仿真分析了基于软件化雷达的脉冲综合孔径雷达SIAR多天线发射系统,基于美国国家仪器NI公司的USRP通用模块
从1894年甲午战争开始,日本帝国主义在历次侵略战争中,对中国人民犯下的滔天罪行罄竹难书。位于旅顺白玉山北麓的"旅顺监狱",就是日本历史罪证的见证者之一。这篇文章,是我在旅
提出了卫星推进系统自锁阀和电磁阀驱动线路的一种优化设计方法。通过数据计算和图形分析阐明新的驱动线路设计的合理性、高安全性和可靠性,避免了常规线路容易出现的干扰和工
本文以台湾运输物流为主要研究对象,通过对运输、物流关联产业发展之探讨,就运输物流发展之外部环境以及台湾运输物流发展现状进行数据搜集,进一步针对台湾内部运输物流和国
陪审制与参审制都是近现代国民直接参加司法审判的制度 ,但两者的具体表现方式及其效果迥然不同。在司法审判民主化改革的理论研究和实践中 ,正确把握两者的区别和优劣 ,对我
该文在分析双闭环V-M调速系统的基础上,根据动态参数工程设计方法,设计了电流环和转速环,确定了时间常数,计算了系统中ASR和ACR两个调节器的相关参数,校验了调节器参数,校核
本研究立足于实习支教受援地区新疆阿克苏地区沙雅县的基础教育发展现状,着重分析反思了当地基础教育改革所取得的一些进步,更进一步地去探讨了改革中所遇到的现实困境并进行