高原地区基于微波辐射计反演大气廓线的神经网络算法研究

以2007~2018年西宁二十里铺气象站探空资料为模拟样本,利用MonoRTM模式模拟中心频率21.985~58.759GHz的35通道亮温,应用BP神经网络对模拟数据进行反复训练,构建最优反演模型,并以2019年探空资料为测试样本,对比分析了不同季节和不同天气条件下BP神经网络与微波辐射计的反演效果。结果表明：晴空条件下,BP神经网络与微波辐射计在温度反演上效果最佳,水汽密度次之,相对湿度最差,其中冬春季BP神经网络反演效果优于微波辐射计,夏秋季反之;有云条件下,BP神经网络温度反演效果在冬、春和夏季均优于微波辐射计,其水汽密度反演效果在四季均较微波辐射计有明显提升,其相对湿度反演效果在冬、春和夏季均较微波辐射计更佳。晴空和有云条件下,BP神经网络在不同季节反演温度、水汽密度和相对湿度的平均绝对误差和标准偏差均小于微波辐射计,尤其是相对湿度的反演精度提升最为明显。晴空条件下,BP神经网络反演温度廓线在春、夏和秋季效果最佳,反演水汽密度廓线在中低层精度较高,反演相对湿度廓线的精度较差,但基本和探空资料趋势一致;有云条件下,BP神经网络反演温度廓线与晴空时基本一致,较微波辐射计精度更高,反演水汽密度和相对湿度廓线在8km以上效果较好。     

直流激电测深二维反演的若干问题研究

目前,激电数据的二维反演主要针对高密度激电数据,而对于常规的激电测深则很少涉及.据此,针对常规激电测深二维反演所涉及到的一些技术难点进行探讨.在正演过程中,根据实测数据实现有限元网格自动剖分,并采用双网格系统,用相邻节点电导率双线性连续进行有限元模拟,并对模拟的误差进行校正.在反演过程中,利用互换原理和Broyden的方法计算偏导数矩阵;采用与奇异值分解算法精度相当的共轭梯度算法求解偏导数矩阵方程.最后通过对模拟数据和实测数据进行反演试算,验证本计算方法是可行的.     

地震反演技术在煤田勘探中的应用探讨

煤炭企业对地质安全保障的期望越来越高,特别是对查清煤层赋存形态、构造、煤层厚度、煤层顶底板条件、岩浆岩侵蚀范围等方面的要求已远超规范。以岩心、测井为约束条件的波阻抗反演,可有效提高目标地质体的预测精度。Jason反演所导出的结果是一宽频带的反射系数序列和宽频带的波阻抗数据,同时加入了低频分量,使反演结果能正确地反映地层变化规律。实例表明,在波阻抗剖面上,不仅可以预测出煤层厚度的变化趋势,而且可以反映煤层夹矸及围岩的变化规律及煤层受岩浆岩侵蚀的范围。     

突发地质灾害气象预警统计模型与应用

通过引入地质灾害潜势度,突发地质灾害气象预警统计模型解决了雨量判据法不能明确表达地质环境条件的问题,在此基础上,提出了潜势度计算、大气降雨变量设计、预警方程优化3方面的改进方法。为提高潜势度计算的准确性,增加了基于卡方检验的地质环境因子独立性判别步骤;对地质环境因子确信程度初始值和权值计算中确信程度变化量初始值设置进行了规定;并将预警模型中累计雨量修正为更为合理的有效累计雨量。在现有的预警方程的基础上提出了基于联合概率分析的预警方程,避免了其与基本地质认识和物理规律存在一定距离、方程系数的意义不明确的局限。以2009年7月2日发布的24 h雨量预报为例,计算结果表明：虽然预报区域面积从101 008 km²减小到27 553 km²,但是地质灾害点落入预报区的比例从44%增加到62%。这说明通过改进的预警模型理论上更加严密,空间准确率有所提高,空报率有所下降。     

透射波波场重建和介质层厚反演研究

本文以P波和SV波波场比值,提出层状介质情况下透射波波场的反向重建算法原理,并据此给出在已知介质速度值时,求介质层厚度的反演方法。文中还以深源远震记录波形资料给出反演算例,证实本文的算法是可行的。     

基于响应面法的高堆石坝瞬变-流变参数反演方法

如何准确确定原级配堆石体的力学参数是高堆石坝建设中亟待解决的一个关键问题,参数反演是解决这一问题的可行方法之一,传统的参数反演方法因需要进行大量的有限元正分析,其计算工作量大,反演效率较低。响应面法可以有效克服以上问题,但已有方法仅针对堆石体的瞬变参数,未考虑流变参数,不能满足堆石坝长期变形预测的需要。综合考虑瞬变和流变参数,通过构造更加合理的响应面函数,提出了基于响应面法的高堆石坝瞬变-流变材料参数反演方法,大大提高了反演的效率和精度。以水布垭面板堆石坝为例,采用该方法对坝体瞬变和流变参数进行了反演分析。反演结果表明,计算值与实测值在数值和变化规律上总体符合较好,反演结果合理可靠且更加高效。     

Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.     

联合GPS与重力资料反演福建地区地壳运动

联合福建省近几年的GPS基本网、基准站观测数据和重力复测数据,反演福建块体的分块运动参数和主要活动断裂运动状态。结果显示,受欧亚板块、菲律宾海板块和太平洋板块相互碰撞的影响,福建地区块体的运动样式大体呈现顺时针方向的运动,造成北东向断裂作左旋走滑运动,而北西向断裂以右旋走滑为主。其中,长乐-诏安断裂带、邵武-河源断裂带和永安-晋江断裂带是块体内最活跃的断裂带,其左旋走滑速率分别为0.7±1.3、0.7±1.7和0.6±1.8mm/a,九龙江下游断裂带具有0.5±1.7mm/a的右旋走滑速率。     

电阻率二维神经网络反演

由于非线性特性地球物理反演一直以来都是一个比较困难的问题. 近十年来,非线性反演方法如人工神经网络、遗传算法在地球物理数据解释中得到越来越多的应用,但目前基本仍限于一维反演问题. 对于二维反问题,反演参数较多,神经网络反演运用较少. 本文利用BP神经网络优化方法,实现了电阻率二维非线性反演. 与传统线性化的迭代反演比较,神经网络反演能够克服传统方法的不足、获得更好的反演结果.     

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.