基于欧拉反褶积法的无人机航磁找矿应用

无人机航磁测量以其高效、安全、低成本等优点逐渐应用到中小面积大比例尺资源勘探领域。本文针对安徽芜湖市三山区22 km²区域构造刻画和找矿前景预测,开展1∶20 000高精度无人机航磁测量工作;推导了欧拉反演求解过程,并分析对比构造指数选取,使用构造指数为0的欧拉反褶积对航磁数据进行反演求解,获得地下磁异常体构造边界解集信息;将其应用到研究区域航磁数据反演,根据反演得到的边界位置和深度信息圈定了18处断裂,并结合异常特征给出了8.7 km²岩浆岩分布区,根据宁芜区域铁矿成矿特点,圈定岩浆分布区及边界为成矿详查区。     

东昆仑胜利铁铜矿区二长花岗岩锆石 Ｕ－Ｐｂ定年、岩石地球化学特征及找矿意义

东昆仑祁漫塔格－都兰成矿带是青海省重要的矽卡岩型铁多金属成矿带，胜利矽卡岩型铁铜矿床位于该成矿带东段，其成矿作用与二长花岗岩密切相关，矿体产于二长花岗岩与大理岩接触带的矽卡岩中。地球化学数据显示，胜利铁铜矿床二长花岗岩ｗ（ＳｉＯ２）为 ７１．８４％ ～７４．５０％，ｗ（Ｋ２Ｏ＋Ｎａ２Ｏ）为 ８．０３％ ～８．６１％，明 显 富 碱；Ｋ２Ｏ／Ｎａ２Ｏ 为 １．２８～１．６２，相 对 富 钾；Ａ／ＣＮＫ为０．９５～０．９８，属于准铝质高钾钙碱性Ｉ型花岗岩。岩石总体富集大离子亲石元素 Ｋ、Ｂａ、Ｒｂ和活泼不相容元 素 Ｔｈ、Ｕ，相对亏损高场强元素 Ｔａ、Ｎｂ、Ｐ、Ｔｉ等。样 品 稀 土 元 素 总 量 ｗ（∑ＲＥＥ）＝１３２．７３×１０－６～１８１．３６×１０－６，（Ｌａ／Ｙｂ）Ｎ 为７．４０～１１．６７，δＥｕ为０．６７～０．７８，稀土元素配分曲线呈右倾型，具有较明显的轻稀土元素富集、重稀土元素亏损的特征，具 Ｅｕ负 异 常。二长花岗岩 ＬＡ－ＩＣＰ－ＭＳ锆石 Ｕ－Ｐｂ加权平均年龄为（２４０．５±２．８）Ｍａ，形成于中三叠世。综合研究表明，二长花岗岩是阿尼玛卿洋盆向北俯冲，致使壳幔岩浆发生混合，后经侵位、结晶分异形成。      

航空重力探测目标发现率研究

航空重力测量相比地面测量可快速地完成区域性的重力数据采集。为在发现异常目标体的同时,高效、安全地完成航空重力测量,作者提出了一种新的目标体发现率的评价方式,指导最优化测量参数设计。通过可靠弱异常的判定标准定义目标体可发现的临界条件,确定目标体可发现概率。通过理论模型分析目标体参数(埋深、水平尺寸等)、测网间距和飞行高度对发现率的影响。绘制目标体发现率关于飞行高度、测线距和测点距的三维图像。     

多旋翼无人机航磁系统误差综合补偿研究

多旋翼无人机磁通门航磁系统以其安全、稳定、高效等特点可广泛应用于中大比例尺矿产资源探测领域,但由于磁通门传感器存在三轴不严格正交,灵敏度、零偏不一致造成的转向差,且系统作业中存在固定场干扰、感应干扰和涡流干扰,需要进行流程繁琐的标定和补偿测试.本文根据实测数据分析出航磁系统机电干扰主要来自机载设备高频干扰,针对机电干扰高频特性设计相应低通滤波器进行误差处理,并基于Tolles-Lawson模型建立仪器转向误差和飞行平台机动误差补偿模型,根据两者结构相似的特点,建立综合补偿模型,该模型简单,实用性强,可通过野外一次补偿测试求取补偿参数进而对工区航磁数据进行补偿处理.最后,将综合补偿研究应用到辽宁省兴城市夹山地区航磁数据,综合补偿后的处理有效去除了航磁数据中的条带状干扰异常,并与地面磁测数据异常形态具有良好的一致性,验证了该研究方法有效性和实用性.     

东昆仑胜利铁铜矿区二长花岗岩锆石U-Pb定年、岩石地球化学特征及找矿意义

东昆仑祁漫塔格-都兰成矿带是青海省重要的矽卡岩型铁多金属成矿带,胜利矽卡岩型铁铜矿床位于该成矿带东段,其成矿作用与二长花岗岩密切相关,矿体产于二长花岗岩与大理岩接触带的矽卡岩中。地球化学数据显示,胜利铁铜矿床二长花岗岩w(SiO_2)为71.84%～74.50%,w(K_2O+Na_2O)为8.03%～8.61%,明显富碱;K2O/Na2O为1.28～1.62,相对富钾;A/CNK为0.95～0.98,属于准铝质高钾钙碱性Ⅰ型花岗岩。岩石总体富集大离子亲石元素K、Ba、Rb和活泼不相容元素Th、U,相对亏损高场强元素Ta、Nb、P、Ti等。样品稀土元素总量w(∑REE)=132.73×10~(-6)～181.36×10~(-6),(La/Yb)_N为7.40～11.67,δEu为0.67～0.78,稀土元素配分曲线呈右倾型,具有较明显的轻稀土元素富集、重稀土元素亏损的特征,具Eu负异常。二长花岗岩LA-ICP-MS锆石U-Pb加权平均年龄为(240.5±2.8)Ma,形成于中三叠世。综合研究表明,二长花岗岩是阿尼玛卿洋盆向北俯冲,致使壳幔岩浆发生混合,后经侵位、结晶分异形成。     

基于三维构造张量的位场边界识别滤波器

地质体构造边界位置的确定是位场数据解释中的一项重要工作,现有很多基于位场梯度张量数据的边界检测滤波器,但存在识别边界位置模糊且无法均衡深浅地质体异常的缺点.本文定义了位场数据的三维构造张量,并提出基于位场构造张量的边界滤波器.为了同时显示不同振幅异常的边界位置,对新定义的滤波器进行归一化处理.在高阶均衡滤波器的计算中需要计算位场的垂向高阶导数,本文引入一种计算的稳定算法,基于拉普拉斯方程利用位场水平导数求解垂向导数,可减小垂向导数计算中产生的误差.将定义的滤波器应用到合成的重磁数据中证明了新方法相比传统的滤波器能更加清晰、准确地圈定边界位置,而且针对同时含有正负异常的地质情况,可避免产生额外的错误边界.最后将新的滤波器应用到实测的重磁数据的解释中,结果显示基于构造张量的滤波器可更准确清晰地划分出断裂的边界位置,发现更多的构造细节.     

基于光滑粒子流体动力学的流动法则对土体滑坡大变形的影响探讨

光滑粒子流体动力学（SPH）是近年来发展起来的一种纯拉格朗日无网格法,并因其在大变形领域内的优势而受到广泛关注。在进行滑坡大变形分析时,流动法则及剪胀角的选取对于边坡失稳后的运动特性有重要的影响。采用Fortran自行编写了基于SPH的边坡稳定性及失稳后大变形分析程序,然后通过2个经典的算例,讨论了不同流动法则及剪胀角的选取对滑坡大变形分析精度的影响。结果表明：（1）剪胀角的选取对土体失稳后的滑动距离有显著影响,随着剪胀角的增大,土体的滑动速度与距离呈明显增大趋势;（2）在关联性流动法则及非关联性流动法则? =1/2? （? 为剪胀角,? 为摩擦角）条件下,土体在大变形过程中会产生过度膨胀,且运动速度与距离要大于实际情况,而采用关联性流动法? =0时,对于非膨胀土可以得到比较令人满意的结果,但对于膨胀土体会使得运动速度和运动距离过小,不符合实际情况。建议在计算膨胀土大变形时,采用非关联性流动法则且适当考虑膨胀性（? =(0.1～0.2)? ）,可以得到较好的结果。     

非对称双环神经网络系统的稳定性和Hopf分岔

环型结构在神经网络中普遍存在,目前对环型神经动力学分岔研究大多数局限于单环情形.值得注意的是,神经网络由成千上万个神经元耦合而成,这些复杂的神经元网络结构不可能只由一个环形结构来准确表述,因此研究具多环拓扑的神经网络模型更具实际意义.本文提出了一种非对称双环神经元网络模型,选择单环的时滞和为分岔参数,分析了双环模型的稳定性和Hopf分岔.最后给出数值仿真对结论进行了验证.     

三疣梭子蟹精子保存研究

于2004年5月到2004年7月对三疣梭子蟹精子进行了常温、低温和超低温保存研究。利用离子载体A23187人工诱导顶体反应以判断精子活力。常温和低温保存实验的结果显示,4℃下无钙人工海水(calcium-free artificial sea water,Ca2 FASW)保存236h的精子有25%以上的成活率。在4℃,25℃两种温度与Ca2 FASW、天然海水(natural sea water,NSW)两种基础液构成的4种组合中,最优组为4℃、Ca2 FASW,最差组为25℃、NSW。将精子浓度与Ca2 FASW对照组及诱导组比较,结果显示,25℃和4℃下,精子浓度与对照组顶体反应率相关性不显著(P=0.068>0.05,P=0.265>0.05),与诱导组顶体反应率极显著相关(P<0.01)。以二甲亚砜(Dimethyl sulfoxide,DMSO)和甘油(Glvcerol,G)作为抗冻剂进行超低温冷冻实验,72h后解冻,包括对照组在内的所有组均有98%以上的精子不经诱导就发生顶体反应。     

基于归一化特征值的重磁数据边界识别方法及其在西藏古堆地热勘探中的应用

基于重磁数据的地质体边界识别是地球物理数据解释的重要内容。随着重磁梯度张量测量技术的不断发展,基于梯度张量数据圈定地质体构造边界位置的方法也在不断发展。本文定义一种基于梯度张量特征值的归一化边界识别滤波器,可同时进行不同深度目标体的边界位置圈定。为验证本文方法的有效性,将建立方法应用于合成的重磁数据中,相比于传统边界识别滤波器,本文滤波器具有更高的水平分辨率和精度,并且在正负叠加异常的试验中,可有效地避免虚假边界的产生。最后,将归一化特征值边界识别滤波器应用于西藏古堆地热田的地质构造解释,结果表明,基于特征值的边界圈定方法能够更准确地显示地质构造的边界位置,发现更多地质体构造细节特征。