三角网格模型剖分方法的研究与应用

数值波场正演在地震学和勘探地球物理学领域正得到越来越广泛的应用,用一个好的方法来建立一个复杂的地质模型则显得尤为关键。三角网格模型相对于层状模型和矩形网格模型在反映地质界面形态、三角网格剖分个数调节及射线追踪速度等方面具有明显优势。利用三角网格剖分的建模方法,在VC2008编程环境下实现了复杂地质模型的描述。通过对复杂模型的建模测试,表明采用三角网格剖分方法,可解决以下模型难题：①对复杂模型界面的描述;②对复杂逆断层的描述;③对封闭块体的描述;④对模型弹性参数进行描述。     

协同克立格法在广西林旺金矿化探数据处理中的应用

协同克立格法,同时兼具化探数据多元性及克立格法表征空间属性的特点。考虑到地质变量两个以上的空间属性,在化探数据处理中运用协同克立格法,可以进一步提高估计精度。运用协同克立格法对广西林旺矿区中金元素成矿进行预测,经相关性分析显示,在矿区中元素的基本组合是Au、Ag、As、Hg,协同克立格插值以Au为主要变量,伴生元素Ag、As、Hg为次要变量。把插值计算后的结果与传统多元统计方法和普通克立格法的计算结果进行比较,结果协同克立格法得到的估值误差较小,预测精度较高,在成矿元素的预测中具有一定程度的优越性。     

3D Zernike矩在三维地形匹配中应用

本文介绍了三维地形匹配工作原理,提出3DZernike矩在三维地形匹配中应用算法。该方法利用三维地形曲面与其3DZernike矩具有一一对应关系,计算实时恢复得到的三维地形与基准图中参考三维地形的3DZernike矩,比较这两个3DZernike矩之间欧式距离,并以欧式距离最小为最佳匹配原则实现三维地形的精确匹配。该方法充分利用3DZernike矩的性质,将三维地形匹配转换成3DZernike矩匹配。该算法对匹配地形无特殊要求,抗噪声能力强,能实现三维地形的精确匹配。     

一种层状大地瞬变电磁响应正演计算的改进方法

频谱法是瞬变电磁法正演计算中常用的计算方法之一,其关键在于内层含贝塞尔函数的积分计算。对于重叠回线装置的瞬变电磁法正演的内层积分,以往采用在积分区间内寻找贝塞尔函数的零点分布,依次在相邻零点之间采用一般的数值积分法在求得各自的积分值后再叠加的方法进行计算。这种方法精度较高,但效率低。利用贝塞尔函数的大宗量渐近特性,提出一种新的计算方法。计算结果表明,该方法计算效率高、方法简单、精度可靠。     

地震动空间相干函数模型中拟合参数的修正

采用SMART-1台阵第5次和第45次地震的水平向分量加速度记录,首先计算了不同间距台站对的地震动空间相干函数;然后讨论了台站距离对相干函数拟合结果的影响,即某一特定距离的相干函数与所有不同距离的相干函数的拟合结果存在明显差异。为减小这一差异,提出了一种对各个不同距离的相干函数拟合参数进行二次回归的方法,并选用Loh相干函数模型进行了验证,最后给出了基于Loh相干函数模型的拟合参数的修正结果。结果表明本文提出的修正方法将大大提高相干函数模型中参数的拟合精度。      

不同地球模型对负荷潮汐改正的影响

负荷潮汐改正与地壳上地幔结构有密切关系,主要原因是负荷潮汐取决于高阶勒夫数。利用较新的1066A地球模型,对负荷勒夫数进行了重新推算,并依据负荷勒夫数组成重力、倾斜、应变及位移格林函数。发现在近区用1066A模型和经典的G-B模型求出的格林函数值差别较大。据此,推估出不同地球模型对负荷潮汐改正的影响,其中远区估算按解析方法进行,近区由数值法进行。     

含气软黏土的不排水抗剪强度计算模型

含气软黏土广泛分布在世界五大洲,大量离散气泡的存在对软黏土的不排水抗剪强度具有强化或损伤作用,然而目前仍缺乏能综合考虑气相强化与气相损伤竞争机制的含气土不排水剪切强度解析表达式。结合临界状态土力学理论,基于笔者提出的含气土屈服函数,建立了v-lnp¢ 空间下临界状态线截距与土体含气特性(u_{w 0},y₀)的数学关系,进而推导了含气软黏土在三轴应力状态下不排水抗剪强度s_u的理论显式表达式。同时,采用所提出的含气软黏土不排水抗剪强度理论计算公式,分别对比了3种典型的含气软黏土在不同含气特性下的不排水抗剪强度理论预测与试验结果,包括马来西亚高岭土、Combwich 黏土和高岭黏土。另外,对比分析了4种不同的含气软黏土不排水抗剪强度理论计算模型的预测误差百分比,证明了提出的基于改进的含气软黏土本构模型的不排水抗剪强度理论计算式的合理性,能够更加准确地同时考虑含气软黏土的气相强化与气相损伤双重效应。     

水库调节地区东江流域非一致性水文极值演变特征、成因及影响

选用8种概率分布函数,以极大似然法估计函数参数,采用AIC、BIC和AICc模型选择准则选出最优分布函数,系统分析变化环境下东江龙川和河源2站的极值流量特征,并对年最大流量变化规律及其影响做有益探讨。结果表明:降雨和水库蓄水工程是年最大流量显著下降的主要原因。龙川和河源站年最大流量LN2混尾分布拟合最好,变化环境后洪水频率最优分布线型基本保持一致,但流量变小造成分布参数改变已导致分布线型高水尾部特性变缓,相应设计流量偏小。用水文情势发生变化前估计的洪水重现期不能很好描述变化后洪水频率特征。     

基于窗口理论的时频分析方法研究

通过对比短时Fourier变换、连续小波变换、S变换、广义S变换等时频分析方法发现,从窗的维度看,它们具有统一的内积形式,差别在于窗函数中σ(f)的取值。通过分析相应方法中σ(f)随信号频率的变化规律,以及随之产生的时窗与时频分辨率问题,解释广义S变换的相对优越性。将广义S变换应用于四类薄互层模型,研究薄互层时频特性与地质特性的对应关系,指导地层旋回方向的预测。     

大别-苏鲁及其邻近地区基于背景 噪声的勒夫波群速度成像

大别-苏鲁造山带是中国大陆东部地区最重要的构造之一. 为了研究该地区的地壳上地幔速度结构,本文收集了国家数字地震台网和中国区域地震台网的山东、河南、安徽、江苏和湖北等省的144个宽频带地震台连续两年(2009年5月—2011年5月)的水平向地震记录(E分量和N分量)数据,首先对台站对之间E-E,E-N,N-N和N-E分量进行互相关,然后分别对这4个互相关分量采用相位权重叠加法进行叠加,最后旋转到横向分量(T-T)获得勒夫波经验格林函数(EGF);用频时分析(FTAN)方法获得4000余条勒夫波群速度频散曲线,并进一步反演得到了周期为6—40s的勒夫波群速度分布图.结果表明,周期为6—10s的勒夫波群速度分布与地表构造特征相吻合.大别造山带、苏鲁造山带、湖北西部隆起均表现为高速;华北盆地发育,表现为大面积的低速;江汉盆地、南襄盆地、合肥盆地等因其规模不同而显示不同程度的低速.在周期为6—30s的勒夫波群速度分布图上,大别和苏鲁地区均显示高速,已有的研究结果中地壳的低速并没有得到反映.其原因一方面可能与勒夫波群速度纵向分辨不高有关,另一方面高压变质岩深度分布可能比已有研究结果给出的要深. 6—30 s的分布图上郯庐断裂带及其邻近地区表现为不同程度的高速,可能与该地区白垩纪以来处于拉张构造体制,地幔物质受到扰动,造成物质上涌有关.