Using constrained inversion of gravity and magnetic field to produce a 3D litho‐prediction model

Geologically constrained inversion of gravity and magnetic field data of the Victoria property (located in Sudbury, Canada) was undertaken in order to update the present three‐dimensional geological model. The initial and reference model was constructed based on geological information from over 950 drillholes to constrain the inversion. In addition, downhole density and magnetic susceptibility measured in six holes were statistically analysed to derive lower and upper bounds on the physical properties attributed to the lithological units in the reference model. Constrained inversion of the ground gravity and the airborne magnetic data collected at the Victoria property were performed using GRAV3D and MAG3D, respectively. A neural network was trained to predict lithological units from the physical properties measured in six holes. Then, the trained network was applied on the three‐dimensional distribution of physical properties derived from the inversion models to produce a three‐dimensional litho‐prediction model. Some of the features evident in the lithological model are remnants of the constraints, where the data did not demand a significant change in the model from the initial constraining model (e.g., the thin pair of diabase dykes). However, some important changes away from the initial model are evident; for example, a larger body was predicted for quartz diorite, which may be related to the prospective offset dykes; a new zone was predicted as sulfide, which may represent potential mineralisation; and a geophysical subcategory of metabasalt was identified with high magnetic susceptibility and high density. The litho‐prediction model agrees with the geological expectation for the three‐dimensional structure at Victoria and is consistent with the geophysical data, which results in a more holistic understanding of the subsurface lithology.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

大地电磁三维正演并行算法研究

大地电磁法三维正演算法计算量大,用传统的串行程序计算相当耗时。而三维正演是逐个频率按顺序计算的,并行性好,适合并行运算。结合MPI自身的优越性,在深入分析大地电磁三维正演串行程序实现流程的基础上,确定了并行计算的思路,实现了三维正演的并行计算。通过三个理论模型对实现的三维正演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维正演并行程序运行结果正确,效率提高明显。此思路可为解决其它地球物理超大计算量问题所借鉴。     

模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

新疆皮山地区地震动衰减关系研究

利用2015年7月3日皮山M_S6.5地震序列的波形资料,采用微遗传算法反演得到震源区的5个介质参数,即应力降Δσ为75.95bars,品质因子Q_0为186.33,η值为0.26,几何衰减分段点R_1为72.18km,R_2为139.70km;并由此计算得到的傅氏谱再配合随机相位谱,进而得到地震动时程,建立了震源区地震动衰减关系。由于缺乏强震记录,本文还收集了2008年10月5日乌恰M_S6.9地震、2008年10月6日阿克陶M_S6.2地震、2012年3月9日洛浦M_S6.0地震的强震数据,将计算结果与经验性衰减关系以及3次强震记录进行比较后发现,皮山地震、洛浦地震的强震数据普遍高于经验性衰减关系,乌恰地震的强震数据与俞言祥等文中的衰减关系吻合较好,在远场低于本文的计算结果。     

GA辅助NLS的GNSS-IR土壤湿度反演方法北大核心CSCD

采用具有阻尼因子的函数模型,使用遗传算法(genetic algorithm,GA)辅助非线性最小二乘(nonlinear least squares,NLS)方法对相位参数进行求解。结果表明:1)相较于标准余弦函数模型,该方法的反演相位与土壤湿度的相关系数有较为明显的提升,反演结果也更加稳定,在5°~15°、5°~20°、5°~25°三个高度角范围内的相关系数均大于0.68,不同高度角之间的相关系数差值小于0.07;2)反演精度有不同程度提高,R 2提高5.72%~76.06%,RMSE减小6.12%~24.24%,MAE减小2.7%~28.3%,将该方案所求相位用于多星线性回归模型后平均RMSE减小10%。     

湖北地区地震震源参数研究

为了研究湖北地区地震震源参数特征,本文基于湖北数字地震台网记录到的S波观测振幅谱,通过遗传算法获得湖北地区的介质品质因子和台站场地响应,并在此基础上计算了湖北地区地震震源参数。结果显示,湖北地区介质品质因子Q值随频率f变化的关系式为Q(f)=501.8·f~(0.309);对于大部分台站获得了与其岩石基底相符的场地响应,没有出现明显放大效应;地震震级与地震矩对数间呈线性关系,即lg M_0=10.06+1.093M_L;震源半径与应力降间呈双对数线性关系,即lgΔσ=10.52-2.01lgr;地震矩与拐角频率间整体上呈反相关;地震矩与应力降间关系不显著。     

2017年6月3日阿拉善左旗5.0级地震震源深度测定

2017年6月3日18时11分内蒙古阿拉善左旗（37.99° N,103.56° E）发生5.0级地震,利用近震深度震相sPL测定法、CAP波形反演法、sPn与Pn震相走时差法,进行震源深度精确测定。研究结果表明,3种方法测得的该地震震源深度基本一致,分布范围为8.5—11 km,表明该地震发生在上地壳。     

Source-independent wave-equation based microseismic source location using traveltime inversion

A new source location method using wave-equation based traveltime inversion is proposed to locate microseismic events accurately. With a sourceindependent strategy, microseismic events can be located independently regardless of the accuracy of the source signature and the origin time. The traveltime-residuals-based misfit function has robust performance when the velocity model is inaccurate. The new Fréchet derivatives of the misfit function with respect to source location are derived directly based on the acoustic wave equation, accounting for the influence of geometrical perturbation and spatial velocity variation. Unlike the mostly used traveltime inversion methods, no traveltime picking or ray tracing is needed.Additionally, the improved scattering-integral method is applied to reduce the computational cost. Numerical tests show the validity of the proposed method.     

四川省城市地震灾害脆弱性综合评价研究

城市化进程的不断推进使得城市的地震灾害脆弱性日益加剧,而城市承灾体的脆弱性受复杂因素影响。本文针对评价指标受主观性影响较大的问题,基于人口、工程、经济及社会4个方面,构建了城市震害综合脆弱性评价指标体系;并构建基于实码加速遗传算法优化投影寻踪（RAGA-PP）的城市震害脆弱性评价模型;最后,对四川省21个市、州进行了震害脆弱性评价。结果表明:巴中、南充等地脆弱性较高,成都、攀枝花等地脆弱性较低;经济因素对城市震害脆弱性影响较大;该评价模型能够克服人为主观性,有效可行。