赣杭构造带中段岩石圈电性结构研究

赣杭构造带（简称赣杭带）位于扬子板块和华夏板块的结合部,是中国重要的铀多金属矿成矿带,发育一系列中酸性火山—侵入杂岩;前人对岩浆喷发和侵入的地球动力学背景仍然存在诸多争议。文章对赣杭带中段（盛源盆地）及邻区收集的大地电磁测深（MT）数据,利用MT相位张量分解技术、二维非线性共轭梯度法（NLCG）等技术手段进行处理和反演获得到了研究区可靠的岩石圈电性结构。结合已有资料分析表明：赣杭带中段浅部火山—侵入杂岩带存在东南方向下倾的高导带和深部软流圈高导区且相连,揭示了该地区中酸性岩浆活动受区域江绍深断裂控制,在晚中生代岩石圈伸展作用下沿着区域深断裂上升侵入和喷发形成火山—侵入杂岩;该岩浆活动处在伸展构造环境与同期Izanagi板块俯冲后撤有关。     

接收函数和面波频散联合反演广西北流MS5.2级地震震源区地壳速度结构

华南陆块属于比较稳定的陆块,地震活动不论是强度还是频度都比西部和华北地区弱很多,2019年在广西北流与广东化州交界处发生的M_S5.2中强地震为探讨华南板内中强地震成因提供了理想场所。本文计算了北流地震震源区布设的7个宽频带流动地震台站、120个短周期流动台站和4个固定台站的远震P波接收函数,结合前人在该区域的面波频散曲线,采用线性阻尼最小二乘反演算法联合反演了研究区地壳三维S波速度结构。得到如下结果:（1）研究区S波速度具有明显横向不均匀性,上下地壳速度在断裂附近存在明显差异性,新丰断裂南侧及东侧上地壳整体S波速度较高;（2）石窝断裂附近存在一条NW向连续性较好的低速带,震中区位置10km深处存在明显局部高速异常,并随着深度减小异常逐渐增加,在浅处其与焦林断裂和新丰断裂附近的高速异常汇聚成一条NNE向连续异常带;（3）结合研究区地质调查结果和余震分布,震中区的局部高速异常可能是NW向石窝断裂南西盘发生右行走滑,并受到新丰断裂高速体的限制和互相作用发生了应力积累,从而引发此次中强地震;（4）地壳厚度泊松比反演结果显示研究区存在NE向明显地壳增厚（26～30 km）,...     

排水沟附近农田潜水蒸发计算的渗流力学方法

潜水蒸发计算方法通常依据潜水水位消退规律建立经验和半经验公式,或是依据土壤—植物—大气连续体建立的SPAC模型。但在设置排水沟控制潜水水位的地段,潜水水位的消退受沟渠排水与潜水蒸发的共同影响,需建立特定水文地质条件下的渗流力学算法。结合实际观测统计制度,将潜水蒸发及灌溉回归等变量离散成阶梯函数;将设置滚水坝的排水沟,概化为水位稳定不变的第一类边界,建立沟渠边界控制下含阶梯函数型源汇项的潜水非稳定渗流模型;利用Laplace变换给出模型解析解,并通过数值法验证解析解的可靠性。依据模型的解,建立利用潜水水位动态实测数据逐日计算潜水蒸发强度的递推公式,讨论沟渠边界的影响;以安徽淮北平原的一个实例,演示方法应用过程与步骤。计算农田潜水蒸发的渗流力学方法,原理上有别于基于潜水位消退规律或SPAC模型的现行算法,且新方法可清晰反映边界条件的影响,计算所需参量较少,易于实际应用。     

基于时间卷积网络的地震波阻抗反演

近些年来,深度学习网络的兴起极大地推动了人工智能技术在地震数据处理、反演以及解译等领域的应用.地震波阻抗反演是石油地震勘探领域的一项关键技术,其反演精度在圈定油气储层构造中起到非常重要的作用.提出了一种基于数据驱动时间卷积网络（temporal convolution network,TCN）模型的地震波阻抗反演方法,旨在无需建立初始反演模型,直接利用工区的少量测井标签数据,以地震振幅数据为输入,将波阻抗反演转化为时间序列建模任务,最终输出地下模型的阻抗信息.采用Marmousi2数据集对基于TCN的波阻抗反演模型进行训练、验证和测试,结果显示,在测试集上该模型预测结果的皮尔逊系数和决定系数分别达到97.92%和95.95%,并对远离训练区域的波阻抗信息预测有着良好的泛化性,且在预测时间和预测精度等方面都要明显优于前人的相关研究工作.上述结果表明,TCN时间序列深度学习模型在复杂地层波阻抗反演中具有一定优越性和应用前景,为地震波阻抗反演提供了新思路.      

新疆西准噶尔苏云河钼矿含矿岩体重力异常特征及其地质意义北大核心CSCD

近年来,新疆西准噶尔巴尔鲁克地区发现的苏云河钼矿、加曼铁列克得铜矿、石屋铜金矿等斑岩型矿床,为研究巴尔喀什斑岩型成矿带是否延伸至我国新疆西准噶尔一带提供了新的研究方向。其中苏云河钼矿是该区最大的斑岩型矿床,赋矿岩体为二长花岗斑岩,目前对于该含矿斑岩体的空间形态及规模认识不足。本文以西准噶尔巴尔鲁克山一带区域重力为研究对象,针对苏云河钼矿含矿岩体开展重力三维反演,从构造环境、含矿岩体、成矿年龄等方面对比分析苏云河钼矿与哈萨克斯坦阿克斗卡铜矿的相似性。结果表明,苏云河含矿斑岩体深部最长约9.4 km,最宽约3 km,垂向延伸超过2.2 km,直观展示了苏云河含矿岩体的规模及空间形态,巨型岩体的存在表明深部仍有广阔的找矿空间。此外,巴尔鲁克山一带的重力场反映了西准噶尔洋残余盆地的特征,证实该区在晚石炭世处于俯冲拉伸相关的岛弧环境,为哈萨克斯坦巴尔喀什成矿带向东延伸到我国新疆西准噶尔一带提供了地球物理证据。     

基于星载激光雷达ICESat-2的雪深监测方法研究

目前,被动微波数据是积雪深度反演的主要数据源,受其较粗空间分辨率的影响,反演雪深存在较大的不确定性。激光雷达由于其较高的测高精度,在雪深监测方面具有一定的潜力,基于星载激光雷达ICESat-2数据对北疆地区2018年10月至2019年9月积雪季的雪深进行了提取。由于很难获取ICESat-2轨迹点的雪深观测资料,因此首先利用地面雪深观测数据对目前流行的被动微波雪深反演产品进行验证,获取精度可靠的雪深产品并与发展的ICESat-2监测雪深数据产品进行对比。结果表明：AMSR2雪深产品在北疆地区误差较大,整体存在高估现象,中国雪深长时间序列雪深产品精度相对可靠,以作为对比ICESat-2模拟雪深的参考数据;ICESat-2雪深与中国雪深长时间序列雪深产品在空间上以及变化趋势方面吻合度较高,但ICESat-2雪深变化更加连续,说明ICESat-2激光雷达数据不但可以提取区域积雪的深度,对积雪深度的空间变化也比被动微波数据更加敏感,可以获取更加详细的积雪深度空间变化细节,为精细化的积雪深度空间分布提供数据支撑。     

基于叠前同时反演的致密砂岩储层预测及含气性识别——以苏里格S区块为例

苏里格S区块储层与围岩阻抗差异小、气水关系复杂,叠后纵波阻抗反演方法难以区分优质储层,通过研究区岩石物理分析可知,叠前参数纵横波比可以有效区分岩性和含气性。本文首先进行模型正演确定储层地震响应特征;然后通过适用于砂泥岩的Xu-White模型进行岩石物理建模及横波预测,建立岩石物理模板;其次采取线性去噪、剩余振幅补偿等处理方法对CRP道集进行优化;最后通过叠前同时反演定量预测储层厚度及含气性。结果表明：(1)储层顶界的地震反射特征为强波谷反射,底界反射不明显;(2)纵横波比小于1.68可有效划分砂岩,联合纵波阻抗小于12 200 g·cm^-3·m·s^-1预测储层,再由更低的纵横波比小于1.57识别含气性;(3)储层展布特征与含气性预测范围趋势相似,但在局部上存在差异,储层发育程度与其含气性不一定呈正相关。本文所提方法期望为下一步圈定有利含气储层面积及井位部署提供有力的技术支撑。     

基于梯度投影法的全变差正则化全波形反演

为降低地震全波形反演的不适定性,常用方法是引入未知模型的先验信息,从而将反演问题正则化。但是,传统正则化方法在包含多个先验信息的情况下,仍然面临挑战。本文提出一种扩展的全波形反演公式,其中包含对模型的凸集约束。本文以慢度平方作为反演的模型参数,展示了如何在施加全变差约束的同时,施加边界约束令其保持在一个物理意义上的可行范围内。为验证本文所提算法的适用性,分别开展简单模型及国际标准地质模型数值实验研究,结果表明,全变差正则化的引入可以提高光滑背景模型下高速扰动体的重构效果。     

大定源瞬变电磁一维自适应正则化反演

自适应正则化反演对初始模型要求较低,直接对瞬变电磁响应数据进行反演。正则化因子通过计算每次迭代的数据目标函数和模型目标函数自适应的得到,从而快速获得地下的地电结构。本文采用自适应正则化反演算法对大定源回线瞬变电磁一维层状模型进行反演,使用均匀半空间作为初始模型,首次采用模型对深度的二阶导数极小的模型约束,通过典型理论模型的反演计算,证明了TEM自适应正则化一维反演算法拟合效果好、精度高,由于不用反复搜索正则化因子,并且收敛速度快,体现了良好的稳定性和可靠性。     

应用随机反演技术精细刻画河道砂体——以大庆油田萨北开发区萨二组为例

以大庆油田萨北开发区萨二组为例,研究了随机地震反演技术及其对河道砂体的精细刻画。研究表明:纵向上,不同厚度砂体的随机反演结果符合率差异较大:厚度大于3 m的砂体符合率为89.6%,厚度为2~3 m的砂体符合率为79.1%,厚度小于2 m的砂体符合率为71.5%;平面上,各沉积单元平均的砂体符合率能达到75%以上,且分流平原亚相的符合率要高于内前缘亚相。地震、测井及地质资料的充分结合,能够对河道的走向、复合河道砂体的单一河道边界及井间砂体叠置关系精细刻画,显示了井震结合预测河道砂体的适用性。