济南泉域岩溶地下水三维等参有限元数值模拟

基于济南泉群研究区域岩溶水贮存和运动特征,建立了三维有限元数学模型,并利用等参有限元方法求解,再根据2001年、2002年的水文资料,识别模型参数,进一步对济南泉域内岩溶地下水的开采方案进行评价,探讨了济南市“保泉”与“供水”的可行性,为科学管理济南市泉群区域内地下水资源提供依据。      

陷落柱的地震识别技术及其应用

通过地震资料解释的研究与实践,对华北石炭二叠系中陷落柱的地震识别技术及其应用情况,进行了较系统的分析和总结。这对今后正确分辨和圈定陷落柱具有一定指导意义。地震剖面、水平时间切片和沿层振幅切片、相干体技术、波阻抗反演技术等,在对陷落柱的发育部位及其空间关系的解释方面,各有其独特效果。山西某矿陷落柱相当发育,利用陷落柱的综合地震识别技术进行解释,效果明显。      

小波分析与C3相干体在小断层解释中的应用

利用小波分析的分频特性及C3相干数据体比常规振幅数据体具有更高的构造分辨能力的特点,将小波分析与C3相干体两种方法相结合,用来识别煤矿采区小尺度构造,特别是识别小断层.实际应用表明,利用小波分析与C3相干体可更加精细地揭示断层等构造信息,提高小断层识别能力.     

形变单位在地震文献中的正确使用

针对地震观测技术论文中有关形变单位"毫秒",即10~(-3″)混用为时间单位ms的情况,笔者从单位符号的意义及观测技术仪器的工作原理出发,最终明确,地震观测仪器中倾斜仪观测量单位应为10~(-3″),不能以时间单位ms代替,二者不能混为一谈。     

基于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方向最为迟钝,最后对以上现象给出了解析.     

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

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

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

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