基于MPI并行算法的电阻率法多种装置数据的三维联合反演

由于地表电阻率法受到浅地表局部异常体的干扰,反演精度受到影响,井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰,但是由于钻井位置的局限性,数据量得不到保障,也会导致反演精度降低.为此,本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先,利用有限差分法实现正演模拟,采用非线性共轭梯度法（NLCG）恢复电阻率结构;其次,调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算,当开辟12个进程时,反演程序获得最大加速比4.51;最后,通过合成数据和实测数据算例证明该反演算法的有效性,也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果.     

THE APPLICATION OF IMAGE-BASED MODELING IN PALEOEARTHQUAKE TRENCH STUDY

The geological structure exposed by paleoearthquake trenches is the key material to the right cognition of fault activity and paleoearthquake. However, paleoearthquake trenching inevitably destroys active tectonic geomorphic evidence and trench exposures are usually difficult to reserve. The conventional process of recording the delicate geological information, manually constructing photomosaics by image-editing software, is time-consuming and produces undesirable artificial distortions. Herein, we explored the process of constructing trench orthophotomosaics and the 3D image model using the Image-based Modeling technology and applied it to the Liutiaohe trench across the Tianqiaogou-Huangyangchuan Fault, Gansu Province. Based on the 3D image modeling and orthophotomosaic, we firstly constructed the control points and scale bars on cleaned trench walls and collected photos of all sections of the trench with a digital camera in the field, and then reconstructed the 3D model of the trench through the Agisoft PhotoScan, an efficient image-based modeling software, and finally yielded the 3D image model of the trench and othophotomasaics of the trench exposures. The results show that the automated workflow can produce seamless, sub-millimeter-level high-resolution photomosaics more quickly, with precision in the centimeter range, and the 3D image model is of great help to identify strata and geological structures in trenches with much lower capital and labor costs and low expertise levels compared with LiDAR, meanwhile, the 3D archive benefits the share and communication and even allows future reinterpreting the site using new insights.     

A NEW FINDING OF SURFACE RUPTURE ZONES ASSOCIATED WITH THE 1936 LINGSHAN M6(3/4) EARTHQUAKE,GUANGXI, CHINA

On April 1, 1936, an M6(3/4) earthquake occurred on the Fangcheng-lingshan Fault. This event is the biggest historical earthquake on the coastal seismic zone, South China ever. But so far, no any findings about the surface rupture of this event have been reported. This paper is the first to find several intact surface rupture zones associated with the 1936 Lingshan seismic event, in the areas of Gaotang, Jiaogengping etc. on the northeast segment of the Fangcheng-Lingshan Fault. According to the field work, the surface rupture stretches to 10km and distributes along NE direction in front of Luoyang Mountain, represented by earthquake scarp, extensional fracture, dextrally faulted gully and river system etc. The characteristics of surface ruptures and faulted landforms indicate that the surface rupture is of normal-dextral strike slip faulting. The trenching on this fault exposed that at least three seismic events have been recorded, including two historical earthquake events and the latest one is the 1936 Lingshan M6(3/4) earthquake. These surface rupture zones are the key to the detection of seismogenic structure and the re-estimate of magnitude of this event. The new finding of these surface rupture zones would be particularly significant for the detection of the seismogenic structure of Lingshan M6(3/4) earthquake.     

CSAMT法三维分辨能力及其物理机制

随着CSAMT法探测应用需求的日益增多,其分辨能力及其物理机制的研究成为人们不得不面对的问题.在CSAMT勘探中,几何效应、电流效应和感应效应是引起电磁场变化的三种物理机制,本文详细介绍了与地下电性结构密切相关的电流效应和感应效应.首先,利用CSAMT法对100 m×100 m×100 m的异常体进行了探测,研究表明存在两个畸变区域.利用电流效应对存在的两个畸变区域进行了解释,进一步研究表明当收发距为10 km时电流效应占据主要地位,感应效应位居次要地位,当收发距增大时电流效应比重下降,感应效应比重增大.选取坐标原点位于电偶极子中心,x方向为电偶极矩方向,z轴垂直向下,随后,研究了CSAMT法对异常体沿x、y及z方向长度变化时的探测灵敏度,研究表明CSAMT法对异常体沿x方向的变化最为灵敏、y方向变化次之、z方向最为迟钝,最后对以上现象给出了解析.     

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.     

基于不等式约束的三维电阻率探测混合反演方法

三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果.     

一种数字三分量地震微测井仪器设计

针对当前地震微测井仪器在施工和数据采集方面的难题,从提高数据采集质量及仪器功能性方面提出了几种改进方法.主要利用MEMS加速度传感器结合多路并行处理采集方法设计了数字检波器,并集成姿态检测技术直接实现水平分量数据的校正功能,保证系统采集的质量;提出了Ethernet-RS485分布式总线自动控制技术,实现各级联探管的随机组合,同时采用电动方式控制探管操作,提高系统的集成化程度.文中详细分析了利用姿态检测技术实现水平分量数据校正的方法,解决了级联探管随机组合控制和震源激发方位的判定问题.仪器单节探管支持三分量并行采集,有效分辨率达18位,参数可远程配置,同时具备快速显示探管姿态的功能,且支持远程控制,非常适于复杂环境下的工程勘探,整个系统实现全数字化操作且仅需四芯电缆互连,从而缩小了设备的体积和重量,降低施工布线的复杂性.     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

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

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

The use of 3D modeling for the prediction of the seismic demands on the gravity dams

Seismic behavior of gravity dams has long been evaluated using a representative two‐dimensional (2D) system. Formulated for the gravity dams built in wide canyons, the assumption is nevertheless utilized extensively for almost all concrete dams due to the established procedures as well as the expected computational costs of a three‐dimensional model. However, a significant number of roller‐compacted concrete dams, characterized as such systems, do not conform to the basic assumptions of these methods by violating the conditions on canyon dimensions and joint‐spacing/details. Based on the premise that the 2D modeling assumption is overstretched for practical purposes in a variety of settings, the purpose of this study is to critically evaluate the use of 2D modeling for the prediction of the seismic demands on these systems. Using a rigorous soil–structure interaction approach, the difference between the two and three‐dimensional response for gravity dams was investigated first in the frequency domain for a range of canyon widths and foundation to dam moduli ratios. Then, the time domain differences between the crest displacements and the maximum principal stress were obtained using 70 different ground motions in order to show the possible bias introduced into the analysis results due to the modeling approach. The results of the study show that even for relatively wide canyons, the 2D analysis can lead to misleading predictions. Copyright © 2017 John Wiley & Sons, Ltd.