基于Tikhonov正则化的双频电磁波电导率成像反演

本文将Tikhonov正则化方法与active-set算法相结合,利用双频电磁波电导率成像原理,求解其反演成像方程.不仅对现有算法进行了改进,也促进了算法的实际应用.本文研究了在双频电磁波电导率成像方程建立后,如何根据其严重病态性质,选择合适的算法求解矩阵成像方程.针对电导率非负的特性,引入正则化参数,将问题转化为一个非负最小二乘问题,并用active-set算法求解.采用改进后的迭代算法对理论模型进行了数值模拟计算,验证了该方法的有效性.应用到实际电导率成像反演,与常规的LSQR、SP-LSQR、Tikhonov正则化等算法进行比较,取得了满意的结果.     

轴对称地层中高分辨率阵列侧向测井信赖域反演法

本文研究轴对称地层中高分辨率阵列侧向测井(HRLA)的多参数信赖域反演方法.首先对前期HRLA的有限元正演方法进行改进,提出基于叠加原理和并行直接稀疏求解器PARDISO的快速正演方案,更适合于反演计算.将HRLA反演问题转化为非线性最小二乘问题,利用信赖域算法求解.为提高反演速度,推导了目标函数对优化参数偏导数的具体计算公式.对典型地层模型,与已有文献中Jacobi预条件共轭梯度法(JCG)计算结果比较,发现PARDIDO比JCG快10倍以上,验证了本文正演程序的正确性与高效性.利用信赖域算法求解了电阻率四参数反演和传统的三参数反演.研究结果表明:并行直接稀疏求解器PARDISO能有效求解此类HRLA正演问题,对6次不同探测深度的测井模拟,所形成的有限元刚度矩阵完全相同,只须进行一次矩阵分解,大大加快了正反演的速度.信赖域算法收敛速度快,且具有全局收敛性.HRLA的信赖域反演结果几乎不依赖于初值的选取,从较差初值出发仍能得到满意的反演结果.另外信赖域算法抗噪能力比较强,即使对测井数据添加信噪比为10dB(甚至5dB)的高斯白噪声,仍能通过反演得到较为准确的地层参数.     

基于非均质边坡强度折减法的三维桥基边坡稳定性分析

基于场变量的强度折减法,以三维一般边坡和陡边坡作为算例模型,对非均质边坡的折减范围进行研究。研究中借鉴某学者研究成果,在ABAQUS中对岩体材料的抗剪强度参数( )和抗拉强度参数T同时进行强度折减,并结合目前边坡失稳判据存在的不足,提出采用坡体内关键部位特征点位移增量 与强度折减系数增量 之比值 和强度折减系数 关系曲线的突变点来确定边坡的临界失稳状态。结果表明,对于一般边坡,只有当黏聚力相差不大的情况下,对黏聚力较小土层实施局部折减才能得到合理的安全系数,对于坡体内存在岩质相差较大的边坡建议采用整体强度折减,相反对于高陡三维边坡而言对相对较软的岩体进行局部折减可以得出较为可信的结果。通过与传统特征点位移法判据相比较,验证了 法的实用性和合理性,同时基于该判据计算出桥基边坡运营期间的安全系数,可为岸坡治理提供理论依据。。     

滇东北峨眉山玄武岩中两阶段自然铜矿化的40Ar/39Ar与U-Th-Pb年龄证据

滇东北与峨眉山大火成岩省有关的自然铜矿化分布于鲁甸、茂林和迤车向斜的二叠纪玄武岩与上覆宣威组地层之间.与自然铜矿化有关的蚀变矿物--浊沸石 (LD- 14- 3、 YSD- 1和 YC- 1)均给出了一致的 40Ar/39Ar坪年龄和等时线年龄 (226～ 228 Ma).这类矿物不存在 Ca和 Cl的干扰,是合适的 40Ar/39Ar定年对象.但它们的 Ar封闭温度较低,在 100～ 110 ℃之间.所得年龄代表了第一次成矿作用的年龄.片沸石的 Ar封闭温度低于 70 ℃,无 Cl干扰的片沸石 LG- 2给出了 (134.0± 1.7) Ma 的 40Ar/39Ar稳定坪年龄, 并与 U- Th- Pb等时线定年结果相一致,表明该区在白垩纪早期存在第二次低温热液作用和自然铜矿化.因此无氯沸石 40Ar/39Ar定年可广泛应用于低温变质作用和冷却热历史研究. 38Ar含量高的阳起石和片沸石样品存在 Cl干扰对 40Ar/39Ar定年的影响.应用 40Ar/38Ar- 39Ar/38Ar等时线关系或 39Ar/38Ar- 36Ar/38Ar相关性进行 36Ar修正有可能获得有意义的年龄值;修正后的阳起石年龄为 235.7～ 238.6 Ma,片沸石年龄为 149.1 Ma,可与无氯沸石的年龄相比较.     

龙门山造山带早期断裂活动的古地震制约——来自汶川科钻(WFSD)岩心的证据

强地震是断裂活动的表现形式,可以诱发地表沉积层序顶部未固结的软沉积物发生变形,形成新的变形层(即震积岩***)。因此,在连续沉积剖面中赋存的多层震积岩应是断裂活动的直接证据。川西前陆盆地中的软沉积物变形记载了龙门山断裂带的活动信息,对认识龙门山造山带演化历史具有重要意义。本文通过"汶川地震断裂带科学钻探"一号孔(WFSD-1)和三号孔(WFSD-3)连续岩心剖面的岩性分析和构造研究,识别出11段不同深度的液化角砾岩层,它们是地震触发成因的软沉积物变形岩层。11个液化角砾岩段厚度从~20m至102m不等,分布在晚三叠世须家河组二-五段。这些液化角砾岩层记录了龙门山前陆盆地形成过程中晚三叠世断裂活动特征及趋势。这些厚度不等的震积岩粗略指示约2~20万年的地震活动长周期(地震幕),以及约4至70万年的间震期(地震幕的间隔时间),反映了龙门山断裂早期脉动式(幕式)活动特征。从不同段液化角砾岩层分布间隔规律来看,地震活跃期间隔(即间震期)越来越短,显示龙门山造山带断裂活动越来越强的趋势。结合前人地表软沉积物变形研究,我们认为龙门山造山带在晚三叠世经历了多期次的正断-逆冲活动的造山作用(至少经历14个地震活跃期),形成龙门山雏形及前陆盆地。     

2008年汶川地震同震滑移特征、最大滑移量及构造意义

2008年汶川地震（M_s8.0）形成了迄今为止空间上分布最为复杂、长度最大的逆冲型同震地表破裂带。沿约275km长的地表破裂带的同震滑移及其最大滑移量的确定,对认识和理解汶川地震地表破裂过程及其变形机制具有重要意义。我们沿地表破裂带进行了详细的滑移特征考察及其同震位移测量,发现沿映秀-北川破裂带分布南北两个滑移峰值区段,南段以深溪沟-虹口破裂段为中心,以逆冲为主伴随右旋走滑运动为特征,最大垂直位移量为6.0～6.7m,北段以北川破裂段为中心,以右旋走滑为主伴随逆冲运动为特征,最大垂直位移量为11～12m,南北两滑移峰值区段所代表的两次地表破裂事件与地震波数据反演结果一致。通过对北川段破裂带的精细地形剖面测量,以及地震前后对比,在北川县曲山镇沙坝村一组获得该破裂段的最大右旋水平位移为12～15m,最大垂直位移为11～12m,这是目前世界上一次地震产生的最大同震垂直位移,最大斜向滑移量为14～17m,为整个汶川地震地表破裂最大滑移量,是汶川地震的宏观震中。北川破裂段高角度的地震断裂、逆冲断裂面的倒转作用以及具最大滑移量的强烈变形作用是北川县城遭受到最强的地表破坏和地质灾害的主要原因。具有走滑量和逆冲量近一致（走滑水平位移/逆冲垂直位移比值为1）的斜向逆冲作用可能是山脉快速隆升的重要机制。     

青藏高原西北部帕米尔东北缘构造地貌与活动构造研究

新生代以来,印度板块与欧亚大陆的碰撞和持续的汇聚在青藏高原西北部的帕米尔地区造成了强烈的陆内变形,形成一系列典型的构造地貌。文章在卫片解译、DEM数据处理的基础上,结合野外地质、地貌观察与测量,对帕米尔东北缘的构造地貌与活动构造特征进行了研究,取得以下认识： 1)在英吉沙地区,通过测量地貌变形面计算出英吉沙背斜隆起高度约为230m,并利用面积平衡法估算出英吉沙背斜的最小构造缩短量约为110m,参考前人的年代学数据计算出英吉沙背斜在中更新世以来的最低隆升速率约为0.23mm/a,最小构造缩短速率约为0.11mm/a; 2)在帕米尔前缘,乌泊尔断裂为一条伴随右旋走滑分量的逆冲断裂,该断裂的右旋走滑作用错断了古近纪地层及流过断裂的河流,通过测量单次地震造成的水系错断量并参考前人研究的该地区大震复发周期约为1000年,估算出该断裂的平均走滑速率为 4.0～6.8mm/a,并推测断裂开始活动的时间大约在 2.2～3.0Ma以前; 3)对喀什地区构造地貌特征的观察与研究表明,明尧勒－喀什背斜和阿图什－踏浪河背斜可能分别为帕米尔东北缘西昆仑山山前冲断带和西南天山山前冲断带的前缘,该地区以西,帕米尔东北缘西昆仑山和西南天山两大构造系统已经发生了碰撞和拼贴。     

各向同性介质长偏移距地震同相轴动校正

传统二阶动校正方法基于较小最大偏移距与目标层深度比和地震波沿直线传播假设,进行长偏移距地震资料处理时,这些假设不再成立.高阶项动校正公式能提高长偏移动校正精度,文中对几种典型的高阶项动校正方法进行了比较,并提出了优化四阶、优化六阶动校正方法.模型计算表明,高阶项动校正方法能取得较常规动校正方法好的动校正结果,但并非阶数越高动校正精度就越高;在纵向速度变化剧烈时,高阶动校正或优化高阶动校正方法一般不能适用于最大偏移距与目标层深度大于3.5的地震反射同相轴,优化四阶和优化六阶动校正公式由于考虑了无穷大偏移距的影响,具有更稳定、更加精确动校正效果,适合于实际的各向同性长偏移距地震资料处理.     

“森林防火地理信息系统”的建设

“９２０２工程”是政府决策ＧＩＳ的典型代表,本文介绍了黑龙江省“９２０２工程”中“森林防火地理信息系统”建设的意义,阐述了系统的设想思想、内容、功能,说明了建立过程、特点,讨论了系统建设中的几个重点问题。     

阿西金矿床流体成矿的元素地球化学界面及X荧光测量识别试验

阿西金矿床元素地球化学界面的形成,是由区域地球化学高背景场和成矿流体作用的结果,其识别标志有岩石组合标志,构造变异标志和蚀变分带标志,同时,对Ｘ荧光测量方法识别元素地球化学界面的标志进行了有效性试验研究,效果良好。